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1.  Excavation  effects  of  underground  nuclear  e^^losions. 


>  During  the  period  1965-170  USSR  conducted  6  underground  excavation  nuclear  explosions. 
Mostly  all  these  explosions  devoted  to  the  investigation  of  excavation  processes  to  use  them 
later  for  peaceful  purposes.  The  explosion  yclds  varied  in  a  wide  range  from  0.2  kton  to  120 
kt.  Explosion  devices  were  buried  into  rocks  at  the  Semipalatinsk  Test  Site  (STS).  General 
data  on  date  and  place  of  events  are  listed  in  the  following  table: 


Name  of 
event 

Date 

Greenwich 
time,  hjn.s 

Coordinates  (degrees, 
minutes,  seconds) 

N  latitude  E  longitude 

Site  name 

BHL  1004 

15.01.65 

6.00.00.8 

49.56.06.0 

79.00.33.7 

BHL1003 

14.10.65 

4.00.00.2 

49.59.26.3 

77.38.08.6 

BHL  101 

18.12.66 

4.58.00.0 

49.55.28.5 

77.44.50.0 

T-1 

21.10.68 

3.52.00.0 

49.43.40.3 

78.29.10.6 

Tellcem 

T-2 

12.11.68 

7.30.00.0 

49.42.44.8 

78.27.40.8 

Telltem 

BHL  125 

04.11.70 

6.02.59.8 

49.55.21.2 

77.45.44.8 

ESSSlSuiil 

Conunent:  BHL  =  borehole 


The  set  of  instrumental  measurements  was  conducted  during  explosion  events.  Particularly, 
seismic  effects,  airblost,  radiation  enviroiunent,  and  dust  cloud  motion  were  monitored.  The 
explosion  dome  motion  was  recorded  by  movie.  For  more  precise  measurements  of  the  dome 
motion,  special  flash  light  markers  were  embedded  at  the  ground  surface  along  the  line 
through  the  ground  zero  point.  The  spacing  of  markers  was  20  to  50  m.  Flash  light  markers 
were  ignited  3  sec  before  the  explosion  and  glowed  about  40  sec. 

Crater  dimensioned  measured  immediately  by  airborne  photo.  As  soon  as  radiation 
environment  allowed,  ground  topographic  measurements  were  conducted.  Main  parameters 
of  6  explosions  craters  are  presented  in  the  table  below: 


Comments:  BHL  =  borehole.  All  crater  dimensions  are  measured  in  respect  to  the  pre-event 


ground.surface.  For  the  T-2  event  the  length  and  width  of  a  channel  are  denoted. 


Sections  below  present  detailed  description  of  rock  properties  at  the  test  site,  as  well  as 
peculiarities  of  the  explosion  process  and  crater  dimensions. 
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1.1  Borehole  1004  event. 

The  explosion  at  the  Borehole  1004  occurred  January  15, 1965, 6h  0.8s  Greenwich  Mean 
Time  (GMT).  The  device  of  140  kton  yield  exploded  at  the  depth  175  m.  The  mm  goal  of 
the  experiment  was  to  test  the  possibility  of  nuclear  explosion  usage  for  industrial  purposes, 
namely  for  large-scale  ground  excavation. 

The  geologic  structure  of  the  event's  vicinity  includes  Low  Carboniferous  rocks ,  covered  by 
relatively  thin  (5  to  7  m)  sheet  of  Triassic  and  Quaternary  deposits.  Low  Carboniferous 
rocks  consist  of  polimictic  and  tuffogenetic  sandstones,  alevrolites,  carbonaceous  shales,  and 
conglomerates.  Rocks  are  shattered  by  tectonic  disruptions,  and  include  dikes  of  porfiritcs. 
Triassic  and  Quaternary  deposits  consists  of  days  and  loam 

The  explosion  device  was  placed  in  shale. 

The  averaged  exploded  rocks  may  be  described  by  following  mechanical  parameters: 

1.  Density  2560kgm-3 

2.  Longitudinal  seismic  wave  velodty  4800  m  s‘  * 

3.  Compressive  strength  430  to  600  bar 

4.  Tensile  strength  70  bar 

5.  Poisson  coeffident  0.34 

6.  Gas  Bearing  Coeffident  (GBC)*)  0.21 

Two  types  of  ground  water  may  be  distinguished  at  the  test  site:  (i)  alluvial  water  at  sand- 
gravd  Quaternary  deposits,  and  (2)  in-fracture  water  at  Low  Carboniferous  rocks. 

Water  bearing  strata  are  hydraulically  connected  and  have  the  united  surface.  At  both 
horizons  ground  water  is  bittered-salted  due  to  chloridc-sulfate-sodium  salt  content. 

Movie  cameras  to  register  the  event  were  installed  along  three  azimuth  at  distances  5, 1 1  and 
1 7  km.  Approximately  40  msec  after  detonation  the  water  began  to  fountain  from  drillholes 
with  the  velodty  of  350  m/s.  At  the  same  time  one  observed  the  dome  formation  with  a 
diameter  of  600  m.  Initial  upward  dome  vdodty  at  the  ground  zero  was  about  100  m/s.  3 
seconds  later  maTimum  velodty  of  geded  rock  was  as  large  as  140  m/s.  The  diameter  of  the 
gecta  column  increased  with  the  velodty  of  10  m/s. 

At  the  end  of  6th  second  the  quick  expanding  condensation  cloud  formed  at  the  top 
of  the  column.  The  condensation  belt  was  observed  around  the  base  of  the  column  during 
the  next  4  seconds.  Approximately  10  seconds  after  detonation  the  ejecta  column  reached  the 
mMimum  hdght  of 950  m  with  the  diameter  of  800  m. 

The  explosion  creates  the  crater  with  the  gecta  zone  around  (Fig.  1.1).  Averaged 
crater  dimensions  are: 

Rim  crest  diameter  520  m 

Apparent  diameter  (at  the  pre-event  surface)  400  to  430  m 
Rim  crest  hdght  above  the  pre-event  surface  20  to  35  m 
Ejecta  zone  outer  radius  (from  GZ)  647  m 

Ejecta  zone  width  (from  rim  crest)  412  m 

Maximum  ejecta  distance  (from  GZ)  900  to  1000  m 

Apparent  crater  volume  below  the  rim  crest  level  10268  m^ 

•)  Gas  Bearing  Coeffident  (GBQ  is  measured  as  the  mass  fraction  of  all  volatiles  removed 
by  heating  of  a  specimen  to  the  temperature  of  1000  to  1  lOOC. 
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Apparent  crater  volume  below 

the  pre-event  surface  6400 

Ejecta  volume  at  the  lim  66 16 

•  Volume  of  the  structural  uplift  at  the  rim  4095  m^ 


Engineer-geological  study  allows  to  distinguish  following  elements  of  the  crater  and 
deformed  zone. 

The  surface  of  the  apparent  crater  (Fig.  1.1)  constructed  of  gravel,  shuttered  gravel 
and  small  fragments  of  test  site  rocks.  Initial  apparent  crater  slope  angles  were; 

North  slope:  30° 

South  slope  55° 

West  slope  35° 

East  slope  40® 


3.5  years  after  the  event  slope  processes  decreased  slope  angles  correspondingly  to  30o,  40o, 
20o,  and  30o  .These  slopes  are  stable.  Due  to  the  downslide  motion,  1 .7E6  m^  sloped  into 
the  crater,  most  of  the  volume  from  the  South  slope.  Downslides  decreased  the  apparent 
crater  depth  at  the  central  part  by  17  m.  The  internal  crater  volume  decreased  at  22%  and 
became  equal  to  5.7  E6  m^in  comparison  with  the  initial  volume  of  6.4  E6  m^- 
Physical  and  mechanical  properties  of  ejecta  deposits  are: 

Spedfle  density  2.5  to  2.85  g  cm-3 

Volume  density  1.61  to  2.23  g  cm-3 

Moisture  content  12  % 

Porosity  45  % 

Coefficient  of  porosity  0.88 

Coefficient  of  water  saturation  0.44 


Size-frequency  distribution  (SFD)  of  fragmented  material  at  South  and  East  crater  slopes 
was  measured  by  deve  analyzes  of  samples  of  shallow  digging,  as  well  as  by  tedmique  of 
oblique  photoplanimetry  (along  the  water  channel  slopes,  and  at  the  surface  of  ejecta 
deposits).  SFD  analysis  show  that  the  main  part  (93.7  %)  of  ejecta  deposits  consists  of 
fragments  below  200  mm.  Fractions  below  200  mm  may  be  separated  into  2  to  20  mm 
fraction  (gravel&pebble  fraction,  Russian  name  "dresva")  which  presents  36.8%,  and  20  to 
200  mm  fraction  (pebble&cobble  fraction,  Russian  name  "tsheben")  which  presents  34% 
Fraction  above  200  mm  ( boulders)  is  disrupted  more  imiformly:  the  coefficient  of  variation 
is  about  17%  in  contrast  to  124  %  for  fractions  below  200  mm. 

Hydrotechnical  investigations  shown  that  the  crater  rim  filtration  properties  allow  to 
use  it  as  a  dam  keeping  the  head  of  10  m  of  water  without  any  additional  sealing  screen. 

The  excavated  rocks  closed  the  river  valley  creating  an  artificial  dam.  To  allow  spring 
melted  water  to  pass  through,  the  special  channel  was  constructed  the  first  spring  after  the 
event.  The  channel  connected  the  crater  with  the  river  valley.  As  the  result,  two  water 
reservoirs  were  created:  the  internal  one  (the  crater  itselQ  and  the  external  one  (in  the  river 
valley  with  the  crater  rim  as  the  dam).  During  the  snow  melting  in  spring  of  the  first  year 
after  the  explosion  10E6  m^  of  water  was  stored  in  the  external  reservoir,  about  7E6  m^  -  in 
the  internal  reservoir.  About  4E6  m^  of  water  infilled  fractured  zones  beneath  the  crater. 
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12  Borehole  1003  event 


The  explosion  BHL  1003  occurred  October  14, 1965, 04h0.2scc  GMT.  The  device  with  the 
'  yidd  of  1.1  kton  detonated  at  the  depth  of  48  m. 

The  test  was  conducted  to  study  the  excavation  process  and  seismic  effects,  including  the 
influence  of  DOB  to  rdative  intensity  of  various  kinds  of  seismic  waves.  Other  goals  were: 

-  to  measure  crater  dimensions  and  the  size  of  residual  deformation  zone; 

-  to  study  residual  rock  transformation  at  various  distances  from  GZ; 

•  to  study  radioactive  pollution  of  rocks  and  ground  water; 

-  to  study  engineering  properties  of  the  ejecta  rim  deposits. 

Local  Geology  of  the  site.  The  geologic  crossection  shows  that  the  site  is  constructed 
of  Upper  Devonian,  Neogenic  and  Quaternary  sediments. 

Upper  Devonian  sediments  arc  presented  by  mostly  ferriferous  thin-layered  alevrolites  and 
argillites  of  a  deep-cherry  color  with  inter-layers  of  gray  sands.  Above  alcvrolite  fat  dense 
days  arc  placed.  These  days  have  greenish-gray  color.  There  thickness  varies  from  0.5  to  27 
m  (22  m  at  the  ground  zero  point).  Quaternary  sediments  are  presented  by  loams  of  3.7  m 
thickness.  The  peculiarity  of  these  sediments  is  a  high  level  of  a  gypsum  content. 


Averaged  properties  of  rocks  around  the  explosion  device 


1.  Density 

2.  Longitudinal  wave  velodty 

3.  Compressive  strength 

4.  Tensile  strength 

5.  Poisson  coeflldent 

6.  Gas  Bearing  Coeflldent  (GBS) 


2690kgm-3 

3650  m/s 
320  to  700  bar 
80  bar 
0.33 
0.054 


are: 


Ground  water.  Across  all  the  event  site  the  water  bearing  horizon  coinddes  with 
Devonian  sediments.  The  device  installation  borehole  was  tamped  by  natural  ground  water 
and  sand.  The  metallic  plug  was  placed  at  the  depth  of  16  m.  Above  the  plug  the  borehole 
was  infilled  with  quartz  sand. 


General  picture.  The  free  surface  motion  at  the  GZ  began  45  msec  after  the 
detonation.  One  can  derive  the  average  compression  wave  speed  of  1 100  m/s.  Before  this 
moment  (approximately  20  msec  after  detonation)  the  light  flash  w^  observed  at  the 
borehole  mouth.  The  flash  looked  like  a  doud  of  radiating  gases.  The  Ught  flash  lasted 
around  60  msec  and  totally  went  off  after  this  moment.  During  this  time  the  light  source 
uplifted  at  5  m  above  the  surface. 

At  the  gatwff  moment  as  the  free  surface  began  to  move,  one  observed  the  water  fountain  not 
far  from  the  borehole  mouth.  The  shape  of  the  fountain  lost  the  visible  contour  very  quickly. 
It  seems  to  result  from  quick  lost  of  a  hydraulic  head.  About  200  msec  the  water  fountain 
was  observed  above  the  growing  dome  surface  as  a  column-like  formation  of  a  constant 
hdght. 

The  initial  upward  ground  motion  at  the  GZ  was  about  42  m/s.  This  value  remamed  a 
constant  during  150  to  200  msec  after  the  ground  motion  began.  Later  (200  to  250  msec  after 
detonation)  one  observed  the  secondary  acceleration  of  the  central  dome  region.  One  may 
connect  it  with  the  gas  acederation  (in  contrast  to  the  free  surface  motion  after  the 
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compressive  wave  reflection).  This  acceleration  continued  about  200  msec.  After  400  to  450 
msec  after  detonation  the  secondary  (gas)  acceleration  stopped.  At  this  moment  the  dome 
summit  uplifted  at  the  height  of  20  m;  the  upward  velocity  was  about  60  m/s.  Later  the  dome 
'  moved  with  a  constant  deceleration  equal  to  normal  gravity  acceleration  till  the  moment  of 
the  dome  disruption. 

At  the  beginning  of  the  gas  acceleration  phase  (210  msec  after  detonation)  dose  to  the 
borehole  mouth  the  new  fountain-like  column  appeared.  Most  probably  this  column-like 
feature  was  the  jet  of  sand  and  water,  infilled  the  borehole,  together  with  construction 
dements  placed  dose  to  the  borehole  mouth.  The  summit  of  the  feature  uplifted  with  the 
velodty  of  70  m/s.  The  feature  was  well  visible  just  to  the  explosion  gas  breakthrough  450 
msec  after  detonation.  The  breakthrough  was  localized  at  the  central  part  of  the  dome.  It 
seems  that  the  borehole  itself  played  as  a  gas-conducting  channel  when  the  tamp  material 
was  thrown  out.  Initial  velodty  of  breakthrough  gases  was  170  m/s.  Initially  the 
breakthrough  looked  like  several  jets  emitted  from  the  base  of  a  column  of  ejected  tamp 
material.  As  the  gas  outflow  velodty  was  much  higher  than  the  ejected  tamp  material,  the 
later  stopped  to  be  visible  very  soon  after  the  gas  breakthrough. 

Breakthrough  gas  jets  very  soon  transformed  into  a  rather  symmetrical  cloud  around  the 
dome  summit.  As  the  dome  uplifted,  this  cloud  uplifted  also,  steadily  expanding. 
Approximately  0.8  sec  after  detonation  one  observed  the  relatively  bright  shine  of  gas 
products  at  the  upper  part  of  the  gas  cloud.  The  glowing  seems  to  connect  with  the  outflow 
of  the  main  cavity  products  to  the  atmosphere. 

The  Rowing  stopped  after  several  tens  of  seconds.  At  the  same  time  the  gas  cloud  growth 
also  stopped.  The  cloud  uplift  velodty  decreased,  and  the  dome  surface  ovemm  it.  The 
cloud  lost  its  shape.  One  saw  separate  jets  from  the  doud  top:  there  were  pieces  of 
constructions  accelerated  during  the  gas  breakthrough.  Some  construction  details  was 
accelerated  to  80  to  100  m/s. 

The  gas  cloud  and  dust  tails  of  flying  construction  fragments  began  to  drift  with  the  wind. 


The  maximum  ground  dome  uplift  was  about  180  to  190  m.  Separated  construction 
fragments  uplifted  at  400  m.  The  general  picture  of  the  dome  growth  is  well  seen  at  the  Fig. 
1 .2,  where  contours  of  the  dome  and  gas  doud  at  different  time  moments  are  plotted.  Table 
1 . 1  presents  the  time  sequence  of  the  explosion  dome  hdght,  and  gas/dust  cloud  top, 

^max  • 
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Table  1.1.  Explosion  dome  height,  and  gas/dust  doud  top,  versus  time  at  the 

BHL  1003  event. 


— # — 

1 

2 

3 

WM 

M 

8 

10 

12 

T.  sec 

■cni 

.98 

ynm 

wm 

warn 

Km 

Kim 

Ull 

■m 

mm 

mm 

■» 

lEiB 

mm 

mtm 

mtm 

HSU 

Km 

Km 

mm 

Bin 

mm 

Kin 

Km 

Km 

Km 

Km 

muiMi 

#  is  the  frame  number  at  the  Fig.  1.2. 

Fig.  1.3  shows  the  time  dependence  of  the  dome  uplift  for  BHL  1004  (I)  and  BHL  1003  (2) 
events.  One  can  see  on  both  curves  the  secondary  acceleration  of  the  epicentral  part  of  the 
dome,  which  seems  to  be  the  result  of  acceleration  due  to  gas  expansion.  For  the  BHL  1003 
event  the  epicentral  uplift  vclodty  increased  from  42  to  60  m/s  during  approximately  150  to 
200  msec.  450  msec  after  detonation  the  dome  uplift  mostly  stopped  and  further  the  dome 
fall  back  with  the  acceleration  about  10  m  s*^  (dashed  line  at  the  Fig.  1.3.) 

Note,  that  secondary  acceleration  at  the  BHL  1004  event  was  recorded  2.35  sec  after 
detonation.  At  this  moment  initial  velocity  of  the  dome  uplift  decreased  from  102  m/s  to  . 80 
m/s.  Then  due  to  secondary  (gas)  acceleration  the  velocity  increased  to  156  m/s. 

The  BHL  1003  explosion  resulted  in  the  crater  with  following  averaged  parameters 
(measured  2  years  after  the  event): 


Crater  rim  diameter  141  m 

Crater  diameter  at  the  pre-event  surface  124  m 

Ejecta  deposits  diameter  28 1  m 

Crater  depth  below  the  rim  crest  29  m 

Crater  depth  below  the  pre-event  surface  20  m 

Rim  crest  height  above  the  pre-event  surface  9  m 


The  crater  and  the  ejecta  zone  are  mapped  at  the  Fig.  1.4.  The  crater  cross-section  along  the 
I-I  line  (see  Fig.1.4)  are  shown  at  the  Fig.  1.5. 

Surface  of  the  apparent  crater  consists  of  clays  and  Devonian  alevrolite  fragments  of  gravel, 
pebble  and  cobble  fractions.  Slopes  of  the  apparent  crater  have  an  inclination  of  40  degrees 
at  the  upper  part,  decreasing  to  8  to  12  degrees  at  the  lower  part.  The  apparent  crater 
volume  below  rim  crest  is  226E3  m^,  the  volume  below  the  pre-event  surface  is  1 1 1E3  m^. 

The  gecta  deposits  consist  of  fragments  and  boulders  of  clay  with  10  to  30  %  of  rubble 
("tsheben")  fraction  of  Devonian  alevrolite.  The  maximum  ejecta  thickness  is  6.5  m,  the 
avenged  width  is  70  m.  Down  to  depth  of  2.5  to  3  m  one  can  see  slag  and  soot  in  clays. 
Physics-mechanical  and  filtration  properties  of  gecta  arc  closed  to  the  initial  ones  for  clays. 
Average  fragment  density  is  2.72  g  cm-3,  average  gecta  dgx>sit  density  is  1 .8  g  cm-3, 
moisture  content  is  13.2  %,  porosity  is  32  to  46  %. 

The  volume  of  gecta  deposits  is  120.5E3  m^. 

For  2  yrars  after  the  event  the  gecta  layer  was  monitored  to  register  the  subsidence.  The 
averaged  subsidence  was  no  more  than  2%. 
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The  zone  of  residual  deformations  (or  the  zone  of  structural  uplift)  is  situated  under  the 
crater  rim,  covered  by  gecta.  The  maximum  structural  uplift  of  6  m  observed  at  the  radius  of 
120  m,  which  corresponds  the  dome  dimension  at  the  moment  of  its  maximum  size.  The 
volume  of  the  structural  uplift  is  60.2E3  m^. 

The  surface  of  the  true  crater  consisted  of  clay  material  at  the  upper  and  middle  parts,  and 
of  pebble&cobble  ("tsheben")  fraction  of  alevrolite  below.  The  content  of  "tsheben"  fraction 
in  clays  increases  from  top  to  bottom.  All  rocks  keep  their  initial  stratigraphy.  The  density  of 
clay  fragments  is  2.72  g  cm-3,  the  average  volume  density  is  1 .92  g  cm  -3.  Moisture  content  is 
24.4  %,  porosity  is  41.7%,  The  internal  friction  angle  is  23  to  27  degrees,  cohesion  is  about 
0.15  to  0.5  bar.  Physical  and  mechanical  properties  of  clay  remain  the  same  as  pre-event  ones 
for  Neogenic  clays.  The  filtration  coefficient  varies  with  depth  from  zero  value  (for  clay 
layer)  to  4  to  500  m/day  at  the  crater  basement. 

The  true  crater  volume  is  235E3  m^  below  the  pre-event  level.  The  volume  of  fallback  rocks 
is  124E3  m^.  The  fallback  thickness  is  43  m  along  the  pre-event  borehole  position. 

The  rock  failure  zone  consists  of  highly  disrupted  alevrolite's  fragment  s  of  pebble&cobble 
("tsheben"),  gravel&pebble  ("dresva")  and  dust  fractions.  The  thickness  of  the  rock  failure 
zone  is  about  1 2  m  below  the  true  crater  floor  along  the  centerline.  The  filtration  coefficient 
is  up  to  300  m/day.  The  volume  of  this  zone  is  13E3  m^. 

Between  the  true  crater  surface  and  the  alevrolite/clay  contact  surface  the  zone  of  spall 
fracturing  is  situated.  The  rode  in  this  zone  is  finely  fractured  alevrolites.  The  levd  of 
fracturing  decreases  with  a  distance  from  the  OZ.  The  rock  fracturing  here  occurred  due  to 
wave  reflection  from  the  boundary  of  different  rocks.  For  this  reason,  in  contrast  to  the 
failure  zone,  here  fine  dust  fraction  and  slag  material  are  absent. 

Thickness  of  damaged  by  spall  zone  at  the  surface  of  the  true  crater  is  6  to  8  m;  maximum 
radius  is  about  80  m. 

The  highly  fractured  zone  presents  by  intensively  fractured  alevrolites.  Fracturing  occurred 
as  along  old  pre-existed  fractures  and  along  new  directions.  The  drill  core  revealed  fragments 
of  10  to  15  cm  height,  rarely  up  to  20  cm.  This  zone  has  a  thickness  of  21  m  along  the  charge 
borehole,  and  goes  up  to  87.5  m  along  the  pre-event  ground  surface.  The  zone  volume  is 
631E3m3 

In  comparison  with  initial  values,  the  strength  of  alevrolites  decreased  by  25  to  35  %, 
Acoustic  impedance  decreased  by  1 1%. 

The  zone  of  blocky  fragmentation  is  presented  by  fractured  alevrolites.  Fracturing  occurred 
along  pre-existed  fractures:  old  fractures  were  renovated  and  became  wider.  The  drill  core 
from  this  zone  presents  column  of  20  cm  height,  rarely  up  to  40  cm  height.  The  thickness  of 
the  zone  along  the  crater  centerline  is  about  20  m,  horizontal  radius  along  the  pre-event 
ground  surface  is  about  115  m. 

The  filtration  coefTident  in  the  zone  of  blocky  fragmentation  is  6  to  7  m/day.  The  zone 
volume  is  1314E3  m^. 
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Hydrologic  monitoring  during  3  years  after  the  event  shown  the  regional  decreasing  of  the 
ground  water  level.  The  radius  of  a  hydraulic  depression  was  above  1000m.  At  the  distance 
of  310  m  from  GZ  the  level  decreased  at  12  m. 

Total  restoration  of  the  ground  water  level  occurred  2.5  years  after  the  serious 

changes  in  the  chemical  composition  of  ground  water  was  observed.  No  Sr^'^  was  observed 
in  ground  water  beyond  the  ejecta  zone. 

Infilling  of  the  apparent  crater  began  3  months  after  the  event.  Maximum  volume  of  the 
water  in  the  crater  was  500  m^.  At  the  end  of  the  3d  year  mineral  content  in  the  water 
achieved  6  g  per  liter  due  to  evaporation  during  summer  seasons. 


1.3.TheBHL  101  explosion 


The  explosion  at  the  borehole  101  occurred  December  18, 1966, 4  hours  58  minutes  GMT. 
The  device  of  the  yield  80  kton  exploded  at  the  depth  227  m. 

The  borehole  101  is  situated  at  the  same  site  as  the  borehole  1003.  Upper  part  of  the 
geological  section  presents  clays  down  to  the  depth  of  40  m,  underlaid  with  sandstone. 
Rocks  around  the  explosion  device  may  be  characterized  with  following  parameters: 


1.  Density  2610kgm-3 

2.  Longitudinal  seismic  wave  velocity  4000  m/s 

3.  Compressive  strength  650  bar 

4.  Tensile  strength  50  bar 

5.  Gas  Bearing  Coefficient  (GBC)  0. 1 28 


When  a  compressive  wave  reached  the  surface,  the  dome  began  to  form  at  the  epicenter 
region.  The  initial  dome  uplift  was  not  recorded  by  movie  cameras  as  the  flash  light  markers 
died  out  at  the  moment  of  compression  wave  impact.  Poor  visibility  prevented  to  get 
measurable  movie  frames  at  the  initial  stages.  For  heights  20  to  50  m  above  the  imtial 
surface  the  dome  uplift  velocity  was  about  48  m/s.  At  the  BHL  101  event  the  uplift  velocity 
was  larger  than  at  the  BHL  1003  event,  despite  the  larger  DOB  value.  It  seems  to  result  from 
the  difference  at  the  local  scaled  geologic  structure  parameters:  at  the  BHL  1003  site  more 
than  a  half  of  the  overburden  rocks  arc  clays  (25  m  of  clays  versus  23  m  of  sandstone);  at  the 
same  time  at  the  BHL  101  site  clay  layer  had  thickness  equal  to  1/6  DOB  only.  Dome 
profiles  are  shown  at  the  Fig.  1.6.  Contours  marked  with  the  time  values  in  sec.  One  can  see 
that  the  distinct  shape  of  the  dome  was  observed  during  approximately  4  sec.  The  dome 
surface  uplifted  at  the  height  of  120  m  and  practically  stopped.  The  dome  basement  radius 
wqs  around  700  m. 

Approximately  4.5  s  after  the  motion  began,  the  light  emission  zones  appeared  at  the  dome 
surface.  These  lights  seems  to  be  connected  with  the  gas  breakthrough  to  the  atmosphere. 
Simultaneously  with  these  lights,  a  lot  of  jets  appeared  at  the  dome.  Jets  united  to  the  joint 
front.  As  the  result,  the  central  part  of  the  dome  was  totally  destroyed.  The  uplift  velocity  of 
the  destroyed  part  increased  from  zero  to  80  to  100  m/s. 

It  is  hard  to  say  from  the  movie  frames,  how  mudt  was  the  mass  of  rock  involved  into  gas 
breakthrough  motion.  It  is  interesting  to  note,  that  lights,  ^ch  were  visible  during  0.7  s. 
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also  uplifted  as  the  dome  evolved.  During  0.4S  s  after  appcannoc,  their  height  increased  at 
50  to  60  m.  If  the  dome  surface  itself  would  be  in  the  same  position,  the  fracture  gas  outlets, 
the  most  probable  place  of  the  light  emission,  would  be  at  the  same  position  too.  So  we  have 
some  ground  to  condude  that  breakthrough  gases  involved  some  dome  material  into 
motion.  An  additional  evidence  is  that  separated  rock  fragments  were  well  visible  after 
separation  from  the  gas  cloud.  The  fall  down  of  these  rocks  generated  an  essential  base 
surge. 

The  explosion  created  the  crater  with  foOowing  parameters: 

1 .  Rim  crest  diameter  404  m 

2.  Crater  diameter  at  the  pre-event  surface  300  m 

3.  Crater  dq)th  below  the  pre-event  surface  1 9  m 

4.  Rim  hdgbt  above  the  pre-event  surface  1 0  m 

5.  Apparent  crater  volume  5.4E5  m^ 

Fig.  1.7  demonstrates  the  BHL  101  crater  cross-section  and  estimated  cavitiy  size.  One  can 
sec  that  the  crater  is  relatively  shallow,  and  has  a  very  rough  floor.  The  central  moimd  is  at 
the  crater  center;  the  top  of  the  central  mound  is  at  the  depth  of  4  m  below  the  pre-event 
surface.  Most  deep  part  of  the  crater  has  a  depth  up  to  19  m.  The  average  crater  depth  may 
be  estimated  as  10  to  15  m. 


1.4.  TheT-1  event 

The  T-1  explosion  was  conducted  October  2, 1968, 3  h  52  m  GMT.  The  device  of  0.2  kton 
yield  detonated  at  the  depth  of  31.4  m.  The  purpose  of  the  experiment  was  to  get  parameters 
of  mechanical  effects  and  radioactive  wasting  for  the  following  experiment  with  the  group 
explosion  of  3  charges  0.2  kton  each. 

The  test  site  constructed  of  Devon  and  Carboniferous  deposits.  The  test  site  situated  at  the 
North-East  wing  of  the  anticline,  constructed  of  sedimentary  Low  Carboniferous  deposits 
(Figs.  1.8  and  1.9).  These  deposits  consist  of  gray  non-uniformly  grained  quartz-feld^ar 
sandstone  with  interlayers  of  dark-brown  argilite  and  yellowish-brown  alevrolite  with  a  high 
content  of  organic  materials.  One  can  Imd  seams  of  carbonaceous  deposits  up  to  3  cm  of 
thickness. 

The  top  24  m  of  rocks  are  highly  weathered;  one  can  easily  crumble  samples.  Sandstone  at 
80  to  85%  consists  of  terrigenic  poor  sorted  material.  The  terrigenic  material  is  presented  by 
fragments  of  effusive  rocks,  quartz  and  plagioclase.  Argilite  and  alevrolite  are  presented  as 
separated  blocks  of  interlayers.  These  interlayers  include  flat  charred  fragments  of  ancient 
plants.  At  the  surface  Low  Carboniferous  deposits  are  covered  by  (Quaternary  clays  and  . 
sands  with  a  thickness  of  3  to  18  m. 


Parameters  of  rock  around  the  explosion  device  may  be  described  as  follows; 


1.  Density 

2.  Longitudinal  seismic  wave  velocity 

3.  Compressive  strength 

4.  Tensile  strength 

5.  Poisson  coeHident 

6.  Gas  Bearing  CoefGcient  (GBC) 


2540  kg  m'3 
4670  m/s 
500  to  700  bar 
140  bar 
0.31 
0.128 
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Underground  water  a  free  surface  at  the  depth  of  2  to  3  m. 

First  motion  at  the  surface  after  the  detonation  was  the  outflow  of  seal  materials  from  the 
charge  borehole  with  the  velocity  of  100  m/s.  Later  the  surface  up  to  radius  of  100  to  110  m 
from  OZ  began  to  uplifr  with  the  initial  vdocity  (at  GZ)  of  43.6  m/s.  Around  the  charge 
borehole  one  observed  the  uplift  of  construction  waste  material  (sand,  gravel).  This  material 
uplifted  5  m  above  the  surface  in  0.4  sec.  The  growing  dome  had  a  distinct  shape. 

The  general  picture  of  the  ejection  process  was  similar  the  event  BHL  1003,  except  no  hot 
gases  breaked  trough  to  the  atmosphere. 

The  evolution  of  the  dome  profile  in  time  is  shown  at  the  Fig.  I.IO.  Table  1.2.  presents 
processed  parameters  of  the  dome  growth  measured  by  photo-  and  movie  registration. 

Table  1.2.  The  dome  growth  at  the  T-1  event.  Time  is  measured  from  the  beginning  of  the 
surface  motion. 


0.12 

0.34 

1.29 

2.24 

3.19 

4.14 

5.09 

6.09 

5.24 

15,0 

66 

105 

135 

157 

162 

157 

Dome 

diameter, 

m 

75 

100 

125 

150 

165 

180 

Uplift 

velocity, 

m/s 

43.6 

44.3 

53.7 

41.0 

31.6 

23.2 

5.3 

-5.3 

Table  1 .2.  shows  that  between  0.34  and  1.29  sec  (after  beginning  of  the  surface  motion)  the 
uplift  velocity  increases  at  10  m/s  relative  to  the  imtial  velocity.  This  is  the  evidence  of  the 
secondary  gas  acceleration  when  the  dome  get  an  additional  velocity. 

Later  the  uplift  velocity  gradually  decreases  to  zero  (at  T=5  sec  approximately)  and  the 
dome  begins  to  descend.  No  breakthrough  of  hot  gases  and  dome  opening  were  observed  at 
this  experiment. 

The  upper  part  of  the  casing  tube  (20  m  length)  with  the  tamp  material  inside  well  overrun 
the  dome  surface  having  maximum  velocity  of  70  m/s.  It  uplifted  at  the  height  of  320  m. 

The  base  surge  began  to  form  3.2  sec  after  the  begiiuiing  of  the  surface  motion  and 
developed  during  40  sec.  The  average  velocity  of  the  radial  spreading  of  the  base  surge  was 
8.5  m/s,  maTimiim  diameter  was  about  700  m,  maximum  height  was  1 16  m. 

The  explosion  created  the  muter  and  the  gecta  2x>ne.  Figs.  1.11  and  1.12  present  the 
topography  of  the  surface  after  the  explosion  and  the  profile  of  the  crater. 

The  apparent  crater  and  gecta  zone  have  dimensions: 


Crater  diameter  at  the  pre-event  surface  60  to  80  m 

Crater  depth  below  the  pre-event  surface  2 1  m 

Crater  rim  crest  diameter  82  to  108  m 

Crater  rim  height  8  to  19  m 

Radius  of  the  qecta  zone  1 10  to  140  m 

Angle  of  the  crater  inner  slope  35  degrees 
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Maximum  flight  distance  of  rock  fragments  was  around  190  m.  Due  to  the  ^^8® 

.  content  the  crater  radius  was  25  %  larger  than  the  predicted  one,  the  crater  d^th  was  40% 
large  than  predicted.  Several  days  after  the  explosion  the  water  was  filled  with  ground  water. 

1.5.  The  T-2  experiment 

The  explosion  T-2  occurred  November  12, 1968, 7  h  30  m  GMT.  Three  explosion  devi^  of 
0.2  kton  each  detonated  at  three  borehole  (#2305,  #2306,  and  #2307)  31 .4  m  below  surface. 
The  distance  between  charges  was  40  m.  The  main  purpose  of  the  exj^mmt  w^  to  obtain 
data  on  mechanical  effects  and  radioactive  wasting  for  future  industrial  using  of 
.  underground  nuclear  explosions.  - 

The  test  site  constructed  of  Devonian  and  Carboniferous  deposits.  The  test  site  is  sit^ted  at 
the  core  of  an  anticline  of  effusive-sedimentary  Upper  Devonian  rocks,  including  porfmte, 
alevrolite,  and  gravelite  (Figs.  1.13  and  1.14).  The  surface  of  Upper  Devonian  rocks  is 
covered  with  Quaternary  clay  and  sand  of  6  m  thickness.  Upper  Devonian  roclw  are 
weathered  to  the  depth  15  to  18  m.  Below  the  weathered  zone  rocks  are  highly  fracture. 
Fractures  contain  calcite. 

Rocks  around  the  device  may  be  characterized  as  follows: 


1.  Density  26 10  kg  m'^ 

2.  Longitudinal  seismic  wave  velocity  4360  m/s 

3.  Compressive  strength  850  to  990  bar 

4.  Tensile  strength  70  bar 

5.  Poisson  coefficient  0*3 

6.  Gas  Bearing  Coefficient  (GBC)  0.04 


Underground  water  has  a  free  surface  at  the  depth  of  2  to  3  m. 

The  first  surface  manifestation  of  the  explosion  was  fountain  of  water  from  exploration 
boreholes  with  the  initial  velocity  of  90  to  1 10  m/s.  The  Section  of  tamp  material  from 
charge  boreholes  and  overrun  motion  of  casing  tubes  were  not  observed.  The  initial  velocity 
of  the  surface  motion  aj  the  GZ  near  the  charge  borehole  #2306  was  43.2  m/s.  The  groimd 
uplifted  as  a  dome  elongated  along  the  line  of  charge.  The  central  part  of  the  dome  uplifted  a 
little  bit  faster.  At  6.3  sec  after  the  motion  beginning  the  dome  reached  its  maximum  size: 
height  of  the  central  part  was  1 12  to  1 15  m,  length  at  the  base  was  205  m,  width  at  the  base 
was  147  m.  The  breakthrou^  of  hot  gases  was  not  observed.  The  dome  growth  is  illustrated 

at  the  Fig.  1.15. 

Along  the  charge  line  special  flash  light  markers  were  installed  at  the  top  of  well  founded 
drill  bars.  A  high-zoom  movie  camera  was  used  to  measure  the  initial  free  surface  velocity 
using  flash  light  markers.  Table  1.3.  summarizes  data  on  initial  velocities. 
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Table  1.3.  Initial  velodtv  of  the  surface  motion  according  to  flash  light  markers  motion. 


#  of  light 
marker 

Measured  uplift 
velocity,  m/s 

The  time 
moment 
of  measurement, 
msec 

Marker  location 

1 

43.2 

50 

At  the  surface,  2.5  m  form  the  central 
charge  borehole  #  2306 

1 

38.4 

150 

The  same 

2 

74.0 

70 

At  the  bar  casted  to  the  case  tube  of  the 
#2306  charge  borehole 

3  1 

33.6 

100 

Between  #2305  and  #2306  boreholes 

4 

20.8 

100 

25  m  from  the  charge  line 

5 

33.2 

100 

Between  #2306  and  #2307  boreholes 

6 

29.6 

100 

15  m  from  the  charge  line 

Main  dynamic  parameters  of  the  event,  measured  by  photo  registration,  are  listed  at  the 
Table  1.4. 


Table  1.4.  Main  dynamic  parameters  of  the  T-3  event. 


^ -  ■■  '  ■ 

Time,  sec 

Parameter: 

0.15 

0.87 

1.94 

3.0 

4.08 

5.15 

Dome  height,  m,  at: 

#2306  GZ 

6.0 

33.4 

70.0 

96.0 

112 

115 

#2305  GZ 

28.0 

54.0 

72.6 

75.0 

75 

#2307  GZ 

• 

29.8 

63.0 

85.0 

95.0 

95 

Dome  uplift  velocity,  m/s,  at: 
#2306  GZ 

38.4 

34.2 

24.3 

15.0 

2.8 

#2305  GZ  1 

- 

24..4 

16.8 

2.8 

0 

#2307  GZ 

- 

- 

31.0 

20.6 

9.4 

0 

Dome  width  at  the  half  height 
level  above  the  day  surface,  m 

90 

100 

115 

140 

165 

Dome  length  at  the  half 
height 

level  above  the  day  surface,  m 

150 

1 

160 

185 

200 

220 

As  in  a  single  explosion  event  T-1,  the  breakthrough  of  hot  gases  was  not  observed.  The 
dome  subsidence  generated  the  base  surge  with  diameter  of  600  to  800  m  and  height  of  40  to 
80  m. 


According  to  the  data  processing  it  seems  that  the  secondary  gas  acceleration  phase  was  not 
as  effective  as  at  the  single  T-1  explosion.  As  a  consequence  the  dome  uplift  velocity  and 
maximum  height  were  lower  that  during  the  single  T-1  explosion.  The  differrace  of  the  gas 
acceleration  phase  intensity  corresponds  to  the  difference  of  volatile  and  moisture  contents 
(expressed  by  the  Gas  Bearing  Coefficient)  at  the  site  of  the  single  explosion  (GBC  =  13  %) 
and  the  multiple  explosion  (4  %). 

Small  difference  in  surface  uplift  velocities  at  ground  zero  points  and  between  them  shows 
the  proper  charge  interference  at  the  chosen  distance  between  charges. 
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The  explosion  of  3  in-line  charges  created  the  trench  crater  and  the  continues  gecta  zone 
(Fig.  1.16).  The  trench  has  a  shape  of  an  elongated  ellipse.  Side  slopes  of  the  crater  between 
charges  are  almost  parallel  to  the  charge  line. 

The  main  parameters  of  the  excavated  trench  are: 

The  trench  length: 

at  the  pre-event  surface  level 
at  the  lim  crest  level 
The  trench  width: 

at  the  pre-event  surface  level 
at  the  rim  crest  level 

The  trench  depth  below  the  pre-event  surface  level  1 6  ra 

Rim  height  above  the  pre-event  surface  level: 


at  trench  sides  9  to  1 

at  trench  ends  7  m 

Continuos  gecta  width: 

at  trench  sides  95  m 

at  trench  ends  65  m 

Maximum  gecta  flight  distance  from  OZ  170  m 

The  trench  volume; 

below  the  pre-event  surface  level  77E3  m^ 

below  the  rim  crest  level  1 87E3  m^ 


1 .6.  Explosion  in  the  borehole  #125 

The  explosion  in  the  BHL  125  occurred  November  4, 1970, 06  h  59.8  s  GMT.  The  device  of 
19  kton  yield  detonated  at  the  depth  1 51 .3  m. 

The  geological  cross  section  includes  porfirite,  alevrolite,  light-green  clay,  sand-gravel 
deposits  and  loam.  These  rocks  may  be  separated  into  two  groups  according  to  their 
physical  and  mechanical  properties.  The  first  group  consists  a  bunch  of  loose  and  cohesive 
materials  with  thickness  of  19  m  at  the  ground  zero.  The  second  group  includes  porfirite  and 
alevrolite. 

Rock  properties  around  the  charge  are  follows: 


1.  Density  2750kgm-3 

2.  Longitudinal  seismic  wave  velocity  4930  m/s 

3.  Compressive  strength  1 1 24  bar 

4.  Tensile  strength  1 10  bar 

5.  Poisson  coefficient  0.25 

6.  Gas  Bearing  Coefficient  (GBC)  0.12 


The  surface  events  were  observed  as  follows.  After  the  compression  wave  arrival  to  the  free 
surface,  the  dome  began  to  uplift  above  the  GZ  region  (Fig.  1.19).  When  the  dome  raised  30 
m  above  the  pre-event  surface,  its  velocity  was  36  m/s.  The  distinct  dome  shape  was  observed 
during  5  sec.  During  this  period  the  dome  uplifted  at  98  m,  the  dome  base  diameter  was  500 


142  m 

160  m 

61  to  69  m 
90  to  108  m 


m.  Approximatdy  3.5  sec  after  the  surface  motion  begmning,  ^  .f 

produL  occurred  at  two  points  of  the  dome  surface.  The  gas 

the  ground.  The  initial  velocity  of  breakthrough  jets  was  35  m/s.  At  5th  and  6th  seconds 
gas  luminosity  was  observed  at  the  vicinity  of  the  breakthrough. 


The  explosion  resulted  in  the  crater  with  following  parameters: 


1.  Crater  diameter  at  the  pre-event  stirfacc 

2.  Crater  rim  crest  diameter 

3.  Crater  depth  below  the  pre-event  surface 

4.  Crater  rim  height 

5.  Apparent  crater  volume 


194  m 
282  m 

17.5  m 

23m 
3E5  m3 


The  crater  cross-section  is  shown  on  Fig.  1.20.  One  can  see  the  relatively  shallow  cratw, 
which  is  very  similar  to  the  crater  of  the  BHL  101  explosion:  the  similar  central  mound 
formed  at  the  crater  center.  The  crater  depth  here  is  about  9.5  ra.  The  rim  is  non- 
symmetrical:  the  rim  height  varies  from  7  to  23  m.  Possibly  this  asymmetry  is  related  to 
oblique  boundary  surface  between  soft  and  hard  rocks. 


the 
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Fig.  1.1.  Geological  cross-section  of  the  crater 
the  BHL  1004  event.  Dimensions  in  meters. 


Ejecta  deposits  outside  of  the  crater  Pre-event  ground  surface 
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Fig.  1.5.  The  crater  and  ejecta  deposits  for  the  BHL  1003  event  (tlie  cross-section  along 
the  1-1  direction  shown  at  the  Fig.  1.4). 


Legend  for  Fig.  1.5. 

1  -  soil  and  loani 

2  -  gray-greenish  and  brown  clay 

3  -  gray-greenish  clay 
4.  -  brown  clay 

5  -  alevrolite 

6  -  diabasic  porfirite 

7  -  clay  with  ale\'rolitc  fragments  (up  to  30  %) 

8  -  crashed  alevrolite 

9  -  zone  of  rock  failure 
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Fig.  1.6.  The  dome  proflle  evolution  at  the  BIIL  101  event.  Time  marks  in  seconds. 


404 
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Fig. 1.9.  Legend  for  the  Fig.  1.8. 


soil  (mostly  sandy  loam) 


gray  sand  of  variable  grain  size,  polimictic 
with  the  admixture  of  10  to  12  %  of  fine  gravel 


Dark-gray  clay  with  violet  hue  (weathered 
crust  of  argillite) 


Brownish-yellow  sandstone,  fine  grained  and 
weat  hered . 


Int er 1 ayer i ng  of  dark-brown  arg i 1 1 i t e  with 
minor  layers  of  yellowish-brown  alevrolite. 
Thickens  of  interlayers  is  of  1.5  to  3  cm. 


Gray  sandstone,  fine  grained  with  thin 
interlayers  of  dark-gray  argillite  and 
gravelite  (thickens  up  to  2  mm). Angle  of 
layering  is  20°  to  25°  to  the  core  axis. 
Rocks  are  slightly  weathered  and  fractured. 
Fracture  surfaces  are  covered  with 
hydroxides  of  Fe  and  Mn. 


Brownish-gray  sandstone  of  middle  grain 
size,  polimictic  with  the  admixture  of 
carbonaceous  fragments  (up  to  4%,  size 
of  1.5  cm,  rarely  up  to  3  cm). 


Dark-gray  argillite,  dense  with  slight 
banding  at  70°  to  the  core  axis.  The 
rock  is  fractured,  with  rusty  fracture 
surfaces.  The  recovered  core  consists 
of  cylinders  of  5  to  15  cm  length.  At 
the  depth  interval  from  32  to  34  m  large 
rock  fragments  was  observed.  A  large 
amount  of  organic  materials  was  found 
at  this  depth  interval. 

Gray  argillite,  of  a  middle  grain  size 
with  interlayers  of  the  fine-grained  one. 
Angle  of  interlayers  is  about  70°  to  the 
core  axis.  One  can  find  infrequently 
thread-like  carbonate  veins. 
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Rim  crest  contour 


Fig.  1.11.  The  topographic  map  of  the  T-1  crater. 
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Fig.  1.14.  Legend  for  Fig.  1.13. 
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soil  ( sandy  loam) 


■  ■/ ;  :  ty 

■  .  -  -4  /  ■ 

•  F  .•  ••  /  .  . 


coarse  debris  with  rubble  fragments, 
polimictic,  with  the  matrix  of  dens 
gray  lime  sandy  loam. 


Motley  clay,  with  yellowish-brown 
and  cherish-red  patches.  Plastic  one. 


i_r  i_r  -L 
lri_r  u~ 
U"  ir  LT 
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Weathered  porfirite,  of  yellowish-brown 
color,  highly  fractured.  Below  the  depth 
of  11.50  m  the  color  is  cherish-red. 

One  see  an  intensive  carbonat i sat  ion . 


Cherish-red  alevrolite,  cohesive, 
with  rare  interlayers  (5  to  10  cm  thick) 
of  coarse-grain  sandstone  (polimictic 
one  with  siliceous  matrix).  It  is  highly 
fractured  rock. 


Gravelite:  gray,  variable  grain  size, 
slightly  rounded,  polimictic,  with 
cohesive  siliceous  matrix.  The  general 
color  is  cherish-red.  At  depth’s 
intervals  26  to  29  m  and  33  to  37  m 
rocks  are  highly  fractured,  angle 
of  fractures  are  inclines  at  15 
to  20°.  Fracture  surfaces  are  covers 
with  calcite. 
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Rim  crest  contour 


The  outer  boundary  of  ejecta  deposits 


Fig.  1.16.  The  topographic  map  of  the  elongated  1-2  crater. 
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2.  Near  surface  nuclear  explosions 


In  the  period  1953  to  1962  at  the  Semipalatinsk  test  site  12  near  surface  bursts  was 
conducted.  Explosion  yields  was  from  0.3  to  450  kton.  The  general  list  of  events  is 
presented  at  the  Table  2.1 . 


Table  2. 1 .  List  of  near  surface  bursts. 


Event 

# 

Date 

q. 

kton 

HOB, 

m 

Target 

rocks*) 

Average 

density 

1 

12.08.53 

450 

33 

1 

9.5 

1700 

2 

05.10.54 

4.0 

0.5 

3 

5  to  6 

1980 

3 

29.07.55 

1.3 

2.1 

1 

20 

1800 

4 

02.08.55 

11.6 

2.1 

1 

20 

1800 

5 

05.08.55 

1.7 

2.1 

1 

20 

1800 

6 

25.03.56 

5.5 

1 

4 

7 

2110 

7 

16.03.56 

14.3 

1 

5 

3 

2380 

8 

09.09.61 

0.5 

0.5 

1 

20 

1850 

9 

14.09.61 

0.3 

0.5 

.  1 

20 

1850 

10 

07.08.62 

10 

0.5 

5 

3 

2290 

11 

25.09.62 

7 

0.5 

2 

- 

1900 

12 

26.11.62 

1.1 

0.5 

1 

22 

1850 

*)  Target  rocks: 

1  -  soft  soil 

2  -  soft  soil  (clayish  deposits) 

3  -  soft  soil  layer  over  quartzite 

4  -  soft  soil  layer  over  highly  fractured  quartzite 

5  -  soft  soil  layer  over  fractured  quartzite 

**)Hrock  the  depth  of  basement  rocks  (or  surface  soft  layer  thickness).  For  events  ## 
2,  6,  7,  and  10,  basement  rocks  are  presented  by  highly  fractured  quartzite,  fractures  are 
filled  with  water. 

In  tables  q  is  the  explosion  yield,  HOB  is  the  height  of  the  device's  center  of  gravity. 

Table  2.2  summarizes  dimensions  of  craters  and  crater's  rims  for  near-surface  bursts. 
Listed  parameters  are: 

Crater  radius  at  the  pre-event  surface; 

Crater  depth  below  at  the  pre-event  surface; 

Crater  volume  at  the  pre-event  surface; 

Rrim  is  the  rim  crest  crater  radius; 

Hrim  is  the  average  rim  crest  height; 

Reed  is  the  continuous  qecta  deposit  outer  radius; 

Rgu  is  the  radius  of  a  maximum  structural  uplift; 

Hsu  is  the  height  of  the  structural  uplift  (averaged  on  azimuths)- 
Meanings  of  these  parameters  are  shown  for  a  typical  crater  on  Fig.  2.2.  Profiles  of 
craters  along  different  azimuths  and  map  schemes  are  shown  on  Figs.  2.2  -  2.  10  . 
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Table  2.2.  Crater  dimensions  of  near  surface  bursts. 


Event 
•  # 

q. 

kton 

HOB, 

m 

Crate 

r 

radius 

1 

m 

Crater 

depth, 

m 

Crater 

volume 

Rrim> 

m 

^rim* 

m 

Rced> 

m 

^su> 

m 

H.SU' 

m 

1 

450 

33 

11.5 

2.4 

500 

-1.6 

45 

2 

4.0 

0.5 

20.5 

7.9 

4940 

24.8 

1.64 

90 

50 

1.26 

3 

1.3 

2.1 

7.25 

3.7 

240 

9 

0.4 

12.8 

12.8 

0.4 

4 

11.6 

2.1 

26 

11 

10700 

31.8 

2.45 

88 

65 

2.26 

5 

1.7 

2.1 

10.75 

5.8 

840 

12.5 

0.85 

18 

18 

0.85 

6 

5.5 

1 

25 

12.5 

11200 

27.5 

2.9 

37 

50 

7 

14.3 

I 

32.3 

17.6 

24150 

41.8 

3.4 

112 

8 

0.5 

0.5 

10 

5.5 

860 

12.5 

1.25 

17.8 

9 

0.3 

0.5 

6.25 

3.5 

215 

712.5 

0.75 

6  to 

13 

10 

10 

22 

9.9 

5050 

28.6 

2.09 

97 

45 

1.2 

11 

7 

0.5 

21.3 

8.35 

5560 

27.7 

1.9 

130 

46.5 

1.45 

12 

1.1 

0.5 

13.5 

7.25 

1410 

16.8 

1.03 

67 
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0.8 

Comments:  Rj-im  is  the  rim  crest  radius,  Hfim  is  the  rim  height,  Rccd  is  the  continuous 


ejecta  deposit  radius,  Rgu 
structural  uplift. 


is  the  structural  uplift  radius,  Hgu  is  the  height  of  the 


The  data  of  geological  cites  for  some  (not  all)  events  arc  listed  below  in  separated  tables, 
marked  with  the  event  number,  "a"  is  the  longitudinal  seismic  wave  velocity.  Strength 
coefficient  (Protodiakonov  Strength  Coefficient)  is  the  value  of  unconfined  compressive 
strength,  measured  in  special  units:  1  unit  =  100  bar.. 


#1 

Thickness,  m 

a,  m/s 

Density, 
kg  m*^ 

Strength 

coefficient 

Heavy  sandy  loam 

0.8 

230 

1700 

0.8  • 

Sand  \vith  15  %  of  gravel 

1.4 

420 

1700 

0.6  • 

Gravel  with  sand  and  pebbles 

1.9 

400 

1900 

0.8 

Loam  with  pebbles 

1.9 

600 

1800 

0.5 

Pebbles  and  cobbles 

3.3 

1000 

1900 

2 

Shale 

6.9 

2000 

2100 

4 

Strong  sandstone 

3.3 

2200 

2200 

8 

Shale 

>10 

2500 

2200 

6 

##3  and  5 

Thickness,  m 

a,  m/s 

Density, 
kg  m’^ 

Strength 

coefficient 

Loam  with  pebbles 

0.5 

- 

1690 

- 

Sandy  loam  with  pebbles 

1.5 

- 

1880 

- 

Gravel 

2.5 

- 

1870 

- 

Loam  with  pebbles 

1.5 

- 

1670 

- 
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#8,  9,  12 


Sandy  loam  ,  sand  with  da’ 


Sandy  loam  ,  sand  with  da’ 


Sandy  loam  ,  sand  with  da’ 


Sandy  loam  ,  sand  with  da’ 


Sandy  loam  ,  sand  with  da’ 


Shale 


Thickness,  m  I  a,  m/s 


Strength 

coeflident 


0.8 


0.95 


1.10 


1.10 


1.22 


1.22 


1.46 


#7,  10 

Thickness,  m 

Loam  with  gravel  and  pebbles 

3 

Fractured  quartzite 

>15 

Strength 

coeffident 


Highly  sanded  day 


Thickness,  m  a,  m/s 


Strength 

coefficient 


Sandy  loam  with  sand 


Sand  with  gravel 


Clay _ 

Water  saturated  quartzite 


Tliickness,  m 

a,  m/s 

1.1 

350 

3.9 

550 

9.5 

550 

2.5 

550 

>5 

500 

Strength 

coefficient 


0.8 


0.8 


1.5 


1.5 


6.0 


I.oam 


Gravel 


Loam 


Sandy  loam  with  pebbles 


Gravel 


uartzitc 


Thickness,  m  la,  m/s 


Strength 

coefficient 


Loam  with  pebbles 


Sajidy  loam  with  pebbles 


Gravel 


Loam  with  pebbles 


Loam  with  pebbles 


Thickness,  m  a,  m/s 


Strength 

coefficient 


Topographic  measurements  of  crater  rims  and  qecta  deposits  was  done  several  years  after 
explosions  for  some  (not  all)  craters.  These  data  are  summarized  below  as  separated 
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tables,  marked  with  the  event  number.  Tables  contain  azimuth  number,  radial  distance 
from  the  ground  zero,  R,  the  measured  height  of  the  structural  uplift  of  the  pre-event 
surface,  h_su,  m,  and  thickness  of  ejecta  deposits,  h_e,  m.  Close  to  the  rim  crest  qecta  was 
absent  at  the  time  of  measurements,  but  it  not  obvious,  were  they  never  deposited  here  or 
been  eroded  in  the  period  of  the  event  and  topo  measurements. 
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Surface  burst  ill  1 


45 


46 


Fig.  2.1.  Tipical  surface  burst  crater. 


R  -  crater  radius  at  the  pre-event  surface; 
h  -  crater  depth  below  at  the  pre-event  surface; 
Rri»  -  rim  crest  crater  radius; 

Hri»  -  average  rim  crest  height; 

Rc«d  -  continuous  ejecta  deposit  outer  radius; 

R«u  -  radius  of  a  maximum  structal  uplift; 

Hsu  -  height  of  the  structal  uplift. 
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Profiles  of  craters  along  different 
for  events  i  .  T  <; 
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of  crater  along  different  azimuths 
for  event  It  4. 


Profiles  of  craters  along  different  azimuths 
for  event  s  lilt  6 ,  11. 


of  crater  along  different  azimuths 


Fig.  2.9.  Isolines  of  ejecta  thickness  for  surface  burst 
events  #  2,  4,  12.  Dashed  lines  -  crater  contur  at  the 

pre-event  ground  level. 


3.  Excavation  explosions  of  HE  charges 

In  1957  an  extensive  series  of  explosion  experiments  was  conducted  at  the  USSR.  The 
purpose  of  experiments  was  to  study  scale  effects  for  increased  depths  of  burst.  These 
experiments  (especially  the  explosion  of  1000  ton  HE  charge)  for  the  first  time  clearly 
demonstrated  the  gravity  effect  at  crater  size. 

The  program  of  experiments  consisted  of  two  sequences:  (1)  explosions  of  charges  with 
mass  10,  80  and  1000  kg  in  loess-like  loam,  and  (2)  explosions  of  charges  with  mass  100 
kg,  1,  10,  and  1000  ton  in  clay.  Explosions  of  100  kg  charges  used  to  compare  excavation 
efficiency  in  two  different  media. 

Crater  and  ejecta  deposits  zone  were  measured  by  topographic  survey.  Measured  data 
used  to  built  the  crater  and  ejecta  zone  map  and  two  cross-sections  along  perpendicular 
directions.  The  later  used  to  calculate  crater  volumes. 


Movie  cameras  with  a  constant  frame  speed  registered  an  explosion  dome  uplift.  The 
velocity  of  the  uplift  was  measured  also  by  contact  gages  installed  at  different  heights 
above  the  initial  ground  surface.  When  the  uplifting  ground  surface  touched  a  gage,  it 
generated  an  electrical  signal,  recorded  by  oscilloscopes. 

Table  3.1  and  3.3  summarize  main  parameters  of  excavated  craters:  Tabic  3.1  contains 
data  for  explosions  in  loess-like  loam.  Tabic  3.2  shows  data  for  explosions  in  clay.  Fig. 
3. 1  -3.6  displays  crater  profiles. 


Loess-like  loam  has  follow's  mechanical  properties: 

1800  to  1900kgm-3 


density 

porosity 

moisture  content 
Clays  properties  are: 
density 
porosity 

moisture  content 


40  to  50  % 

10  to  30% 

1800  to  2200  kgm-3 
5  to  10% 

25% 


Soil  density  increased  with  depth.  At  the  same  time  near  to  surface  soil  density  was 
sometime  less  than  listed  above:  surface  layers  of  loess  had  density  of  1600  kg  m'^- 

Experimental  study  allowed  to  separate  schematically  the  excavation  process  into  3 
stages. 


The  first  stage  endures  from  the  moment  of  detonation  to  the  moment  when  a 
compressive  wave  reaches  the  free  ground  surface.  The  ground  motion  here  is  exactly  the 
same  as  for  totally  contained  explosion.  The  free  surface  does  not  influence  the  motion. 
The  time  duration  of  the  synunetrical  motion  is  relatively  short,  but  this  stage  is  very 
important  for  consequent  effects.  For  example,  experiments  show  that  explosions  with 
the  characteristic  ratio  n  <  3  the  explosion  cavity  has  enough  time  to  grow  up  to  the  size 
close  to  camouflet  cavity.  This  fact  demonstrates  that  at  excavation  explosions  a  large 
part  of  energy  is  spent  to  plastic  deformations  and  consequent  heating  of  the  target 
material. 
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At  the  second  stage  energy  of  detonation  products  is  used  mostly  to  accelerate  the  ground 
in  the  direction  to  the  free  surface.  The  free  surface  distorted  previously  central- 
symmetrical  velocity  field  and  maximum  uplift  velocity  may  be  essentially  larger  than  the 
particle  velocity  at  the  compression  wave  front.  As  the  result,  the  specific  dome  of 
uplifted  ground  is  formed.  Thickness  of  the  dome  shell  is  small  and  continues  to  decrease 
as  the  dome  uplifts.  Later  the  dome  disrupts  into  a  swarm  of  flying  fragments.  At  the  end 
of  the  second  stage  the  energy  of  detonation  products  is  totally  consumed. 

Third  stage  is  the  inertial  flight  and  deposition  of  fragments  in  the  gravity  field. 

The  1957  experiments  created  the  foundation  to  construct  the  above  mentioned  scenario 
of  crater  excavation  explosions. 


Table  3.1.  Explosions  in  loe.ss-like  loam 


Charge 

mass, 

ton 

Don, 

m 

Crater  dimensions 

n 

Radius, 

m 

Depth,  m 

Volume, 

m 

1 

0.01 

0.6 

1.5 

0.9 

_ ^ _ 

2.5 

2 

0.01 

0.65 

1.7 

0.95 

a) 

2.65 

3 

0.01 

0.95 

1.75 

0.8 

a) 

1.84 

4 

0.01 

1.0 

1.75 

0.8 

a) 

1.75 

5 

0.01 

1.2 

1.8 

0.75 

._il _ 

1.5 

6 

0.01 

1.55 

1.8 

0.7 

a) 

7 

7 

0.08 

2.0 

3.25 

1.8 

27 

1.62 

8 

0.08 

2.5 

3.85 

1.9 

36 

1.54 

9 

0.08 

2.75 

3.7 

1.7 

32 

1.35 

10 

0.08 

3.4 

3.75 

0.5 

17 

1.1 

11 

1.0 

4.0 

8.0 

3.4 

385 

2.0 

12 

1.0 

5.0 

7.0 

2.85 

225 

1.4 

13 

1.0 

6.0 

8.4 

3.95 

425 

i.4 

14 

1.0 

7.0 

8.5 

3.1 

315 

1.2 

a)  Volume  is  negligible 
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Table  3.2.  Explosions  in  clay 
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4.  Physical  Modeling  of  Excavating  Explosions 

The  thorough  study  of  cratering  explosions,  presented  above,  namely  the  scheme  of 
consequent  stages  of  energy  transfer  into  excavated  rocks,  resulted  in  the  special 
technique  of  large  scale  excavation  explosion  modeling. 

Suppose  that  the  first,  shock-wave  stage  of  explosion  has  been  terminated.  The 
dimensions  of  shot  cavity  arc  thus  known.  The  rock  around  the  shot  cavity  has  suffered 
disintegration  free-surfacewards.  The  breakdown  zone  radius  is  also  known;  as  a  rule,  it 
is  2-3  times  the  charge  burial  depth.  For  this  reason,  the  original  strength  properties  of 
the  rock  are  of  minor  importance.  The  shattered  rock  is  visualized  as  a  medium 
composed  of  grained  particles.  The  cohesion  between  separate  fragments  of  this  medium 
is  not  great,  and  the  medium  offers  resistance  to  shear  stress  in  accordance  with  the  law 
of  dry  friction.  No  kinetic  energy  imparted  to  the  rock  at  the  shock- wave  stage  is  taken 
into  account.  Determinative  of  the  process  is  motion  of  broken  material  in  the 
gravitational  field  by  the  action  of  gas  acceleration.  The  work  done  by  gaseous  explosion 
products  is  dependeirt  on  thermodynamic  parameters  and  burial  depth  of  explosion. 

As  we  discuss  above,  the  excavation  process  may  be  described  as  qection  of  crushed 
rock  by  the  gaseous  explosion  products.  The  initial  explosion  stage  within  which  the  shot 
cavity  is  formed  and  the  rock  undergoes  destruction  is  ignored,  and  its  involvement  is 
reflected  only  in  the  choice  of  appropriate  parameters  of  the  initial  model  state.  The 
zero-of-time  for  the  initial  model  state  is  taken  the  time  of  termination  of  the  first  stage. 
Therefore  the  shot  cavity  radius  r,  the  cavity  gas  energy  E  (or  pressure  P),  and  the 
explosion  depth  W  are  taken  as  the  initial  model  parameters.  Since  the  rock  ejection 
process  is  a  slow  one,  the  medium  compressibility  may  be  ignored.  To  characterize  the 
broken  rock,  we  make  use  of  the  density  p ,  the  internal  friction  coefficient  K,  and  the 
cohesion  C  (to  account  for  the  cohesion  of  qected  rock  with  surrounding  massif. 
Naturally,  the  gravitational  acceleration  g  and  the  atmospheric  pressure  Po  over  the  free 
surface  are  to  be  included  among  the  determining  parameters.  The  dependence  of  the 
major  crater  dimensions  (for  example,  radius  iJ)  on  the  initial  conditions  and  the 
properties  of  ejected  rock  reads: 

mv  -  F(EJ p gW*;  E/ Po  ;  E/c  ;  W/r;  k)  (4.1) 

Two  limited  cases  arc  distinguished  in  relationsliip  (4.1).  First,  when  parameter  E/ pgW* 
may  be  left  out  of  consideration,  a  geometrical  similarity  takes  place,  and  the  cavity  gas 
energy  is  E~  . 

.Apparently,  a  necessary  and  sufficient  condition  of  this  is  the  smallness  of  quantity  pgfP 
when  contrasted  Po  and  C.  This  condition  defines  small-scaled  explosions.  In  the  other 
limiting  case,  when  the  scale  of  an  explosion  is  large,  the  only  parametric  term 
incorporating  energy  isE/pglP^ .  ITiis  signifies  that  similar  craters  with  E/hF  =  ermsf 
arc  generated  at  E/ pgW^  =  const,  that  is,  the  cavity  gas  energy  must  be  proportional  to 
,  rather  than  to  as  one  should  expect  on  the  geometrical  scaling  basis.  If  the 
medium  characteristics  remain  invariant,  the  role  of  parameter  E/  pgW*  becomes 
increasingly  more  decisive  with  higher  scaling  of  the  event. 

As  suggested  by  relationship  (4.1),  it  will  suffice,  to  comply  with  the  modeling  scaling 
conditions,  to  diminish  the  dimensional  parameters  Po ,  c,  E,  W,  r  such  that  the 
dimcntionless  combinations  in  the  modeling  experiment  w'crc  the  same  in  value  as  their 
real-life  counterparts.  To  this  effect  with  a  view  to  reduce  the  value  of  /*o ,  an  evacuation 
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technique  has  been  applied  to  the  atmosphere  over  the  free  surface.  In  order  to  reduce 
the  value  of  c,  granular  and  unconsolidated  material  have  been  used  in  modeling,  with 
.  insignificant  cohesion  and  grains  much  smaller  in  size  than  the  explosion  depth  IK  This 
choice  of  the  model  material  was  in  conformity  with  the  schematization  conditions 
which  require  the  ejected  rock  be  unconsolidated.  The  model  cavity  gas  energy  E  was 
brought  in  correspondence  with  the  explosion  depth  W  to  obey  the  scaling  conditions. 
The  model  cavity  radius  r  was  chosen  subject  to  the  condition  of  maintaining  constant 
fV/r  ratio.  Thus,  the  underlying  idea  of  the  adopted  modeling  method  is  transformation 
of  force  fields  aimed  at  maintaining  the  same  proportion  of  forces  operative  both  in 
model  and  in  real  life.  A  controlled  variation  of  the  free  surface  pressure  and  of  the 
material  cohesiveness  will  allow  to  reach  the  needed  relation  of  determining  parameters 
subject  to  scaling  conditional  in  conduction  a  properly  scale  explosion. 

In  particular,  in  modeling  large-.scaled  excavating  explosions  -  a  major  pursuit  of  our 
endeavors  -  the  effect  of  parameters  E/  Po  and  E/c  was  reduced  to  a  linear  case 
owing  to  the  use  of  high  vacuum  and  dry  quarts  sand.  Under  such  conditions,  the 
dimensionless  parameter  E/pglV^  becomes  a  determining  one.  To  meet  these 
conditions,  a  mas.sif  of  dry  .sand  was  loaded  in  a  vacuum  chamber.  Within  the  sand 
massif,  a  volume  of  eompressed  air  encased  in  a  spherical  elastic  envelope  was  emplaced 
at  a  definite  depth.  The  sudden  expansion  of  ga.s  from  rupturctl  envelope  produced  an 
excavatory  action  on  the  overlying  layer  of  sand  to  form  a  crater.  In  this  experiment,  the 
dimensionless  combinations  of  determining  parameters  were  in  correspondence  with  the 
conditions  for  conducting  a  large-scaled  explosion. 

MODELING  FACILITIES.  In  designing  our  physical  model,  a  number  of  versions  have 
been  tried.  Schematically,  a  diagram  of  one  of  this  installation  is  .shown  in  Fig.  4.1.  The 
vacuum  chamber,  shaped  as  a  80*80*100  cm  rectangular  parallelepiped  with  stiffening 
ribs,  was  made  from  sheet  steel  8  mm  thick.  The  front  side  of  the  chamber  was  made 
from  a  transparent  Plexiglas  plate  5  cm  thick  to  allow  \'isual  inspection.  The  upper  part 
of  the  back  side  could  be  hinged  out  to  admit  access  for  the  preparation  of  an 
experiment  and  for  making  post-shot  mcasuremaits  of  the  crater  formed.  Model 
material  was  loaded  in  the  chamber  to  fill  up  about  a  half  of  the  available  volume.  The 
prc.ssure  inside  the  chamber  could  be  set  at  any  level  from  the  atmospheric  down  to  10'^- 
10'<^  kg/cm^.  A  standard  vacuum  gauge  was  used  to  measure  the  subaimospheric 
pressure.  The  "explosive  charge"  was  a  thin  rubber  envelope  of  spherical  shape 
contaming  a  working  volume  of  air;  the  envelope  was  enclosed  within  a  wire  mesh 
framework  which  thus  specified  the  cavity  volume.  The  clear  area  of  the  wire  mesh 
framework  surface  was  80-90%  and  practically  did  not  affect  the  accuracy  of 
experimental  measurements.  During  the  chamber  evacuation,  the  cxce.ssivc  pressure  in 
the  working  volume  was  controlled  at  a  given  level  with  the  aid  of  a  diflerential 
manometer.  The  differential  manometer  was  an  air-tight  cylinder  5  cm  in  diameter  anti 
60  cm  in  height,  filled  in  part  with  mercury.  A  tube  with  an  open  end  8  mm  in  diameter 
ran  dowm  the  mercury  column  nearly  touching  the  cylinder  bottom.  The  upper  end  of 
the  lube  was  connected  to  the  working  volume,  the  space  over  the  mercury  surface 
communicated  with  the  vacuum  chamber.  Making  the  mercury  level  vary,  one  can 
control  a  specified  excessive  pressure  inside  the  working  volume.  In  the  course  of 
evacuation,  the  excess  air  from  the  working  volume  passes  through  the  mercury  to 
mauitain  an  excessive  pressure  in  the  working  volume.  The  working  volume  pressure 
with  respect  to  the  pressure  over  the  material  surface  is  controlled  by  the  mercury 
manometer  and  at  need,  with  the  differential  manometer  shut  off  the  system,  can  be 
raised  to  the  aimo.spheric.  The  needed  vacuum  having  been  reached,  the  envelope  is 
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ruptured  through  its  contact  with  a  Ni-Cr  alloy  filament  heated  with  electric  current. 

The  gas  volume,  suddenly  released,  ejected  the  model  material  to  form  a  crater.  The  time 
of  envelope  rupture  was  automatically  synchronized  to  the  start-up  of  a  high-speed 
instrumentation  camera  which  recorded  the  throwout  process. 

This  installation  has  been  used  to  model  nuclear  cratering  e.xplosions  and  certain 
famiKar  large-scaled  chemical  explosions  under  the  conditions  of  a  simple  free-surface 
topography.  The  experimental  model  results  compared  with  the  results  of  field  test 
explosions  in  both  the  crater  dimensions  and  kinetic  parameters. 

Table  4.1  summarize  values  of  crater  radius-to-DOB  ratios  for  several  large  scale 
explosions  and  their  model  analogs. 


Table  4.1 


#  . 

Explosion  event 

Real  value 
ofR/W 

Model  value 
ofRfW 

1 

BHL  1004 

1.18 

1.2 

2 

BHL  1003 

1.7 

1.12 

3 

T-1 

1.19 

1.15 

4 

Sedan 

0.95 

0.9 

5 

Schooner 

1.2 

1.27 

6 

Cabriolet 

1.06 

1.03 

The  further  elaboration  of  the  modeling  tcehnique  allowed  to  simulate  large-scale 
explosions  with  complex  mountain  environments.  Beginning  1972  all  projects  of 
hydropower  construction  u.sing  HE  explosions  were  pre-tested  by  the  model  technique, 
described  above.  3D  models  of  real  explosion  sites  were  created  in  the  vacuum  chamber 
using  sand  mixture  with  up  to  1%  of  glycerol.  The  power  source  for  the  explosion 
simulation  was  a  pressurized  rubber  bag.  The  modeling  allowed  to  e.stimate  the  optimum 
consequence  of  charges  position,  energy,  and  detonation  sequence  to  create  the  project 
dam  profile. 

The  typical  example  is  the  explosion  construction  of  the  dam  at  the  river  Uch-Tcrck. 

The  explosion  site  was  situated  at  the  narrow  valley  with  slope  angle  up  to  50  degrees 
and  dq>th  up  to  200  m.  Explosion  charges  were  buried  under  the  right  slope  at  two 
levels.  The  low  level  charge  of  544  ton  had  a  length  of  145  m  and  DOB  (normally  to  the 
free  surface)  50  to  65  m.  The  upper  charge  of  620  ton  had  the  length  of  92  m  and  DOB 
of  82  to  94  m.  Charge  masses  are  cited  in  equivalent  tons  of  TOT  as  the  lower  charge 
was  constructed  of  ANFO  (5.5%  of  fuel  oil).  The  upper  charge  was  designed  for  rock 
fragmentation  at  the  upper  part  of  the  valley  slope.  The  lower  charge  detonated  5 
seconds  later  than  the  upper  charge.  It  was  designed  to  push  out  rocks,  sliding  down 
after  the  first  charge  detonated. 

Fig.  4.2  shows  profiles  of  the  constructed  dam  across  and  along  the  valley.  The  dashed 
line  corresponds  to  the  model  predictions.  One  can  sec  the  reasonably  goo  fit  of  the 
model  and  the  real  dam  profiles. 
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Fig.  4.1.  The  modelling  installation  diagram:  1  -  vacuum  chamber; 

2  -  working  volume;  3  -  vacuum  pump;  4  -  vacuum  gauge;  5  -  oil 
pressure  gauge;  6  “  differential  manometer;  7  —  mercury  manometer; 
compound  pressure  and  vacuum  gauge;  9  “  compressed  air  cylinder; 

10  -  blasting  unit;  11  -  instrumentation  camera. 
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Across-  and  along-valley  crossections  of  the 
dam  at  the  river  Uch-Terok 

1  -  the  real  explosion  dam,  2  -  the  model  simulation 


5.  Z-model  for  near  siarlace  bursts. 

5.1.  Historical  notes. 

In  the  early  70th  the  model  oi  the  excavation  ground  was 
created,  named  by  US  researchers  as  Z-model.  Independently,  in 
1973-15  the  similar  model  was  elaborated  in  Russia.  First 
unclassified  publications  seems  to  appear  simultaneously  in  1976 
in  Abstracts  for  the  Symposium  on  Planetary  Cratering  Mechanics, 
Flagstaff,  Arizona,  13-17  September  1976  and  next  year  [4,  5,  6, 
7].  Only  in  1994  Russian  researches  got  a  possibility  to  look 
through  the  one  of  basic  American  reports  described  Z-model  [Sl¬ 
it  become  obvious  that  in  1974  Maxwell  and  Seifert  always  done 
many  of  applications,  which  was  re-invented  in  Russia  during  the 
late  70th.  The  first  unclassified  detailed  description  for  the 
Russian  version  of  Z— model  was  published  only  in  1979  [9]-  So  no 
questions  about  priority  arise-  Nevertheless  is  seems  to  be 
instructive  to  describe  the  Russian  approach  and  philosophy  for 
the  Z-model,  as  well  as  some  Ideas  which  we  never  met  in  American 
unclassified  publications • 

Initially  we  named  our  model  as  Integrally  Balanced  Flow 
(IBF)  [6]  or  Simple  Model  of  Cratering  (SMC)  [7,  9]-  As  we  mostly 
used  the  letter  n  for  the  velocity  exponent,  it  might  be  called 
"n-model”,  but  here  we  will  name  the  model  with  its  real  name. 

5.2.  Main  ideas. 

As  we  can  conclude  from  [8]  the  starting  point  for  Maxwell 
was  computer  simulations  of  surface  bursts,  revealed  the  specific 
geometry  of  the  velocity  field  at  late  stages  of  cratering.  In 
Russia  the  starting  point  based  on  experiments  with  a 
quarter-space  geometry  [9,  10]. 

HE  charges  of  PETN  detonated  at  the  surface  of  wet  and  dry 
sand  in  the  box-  One  wall  of  the  box  was  a  transparent  one 
(plexiglas).  The  steel  plug  of  2  cm  diameter  closed  the  hole  in 
the  Plexiglas  window  to  keep  the  window  from  damage.  Charges 
detonated  at  the  plug  surface.  In  some  of  experiments  the  plug 
was  locked  with  a  narrow  vertical  steel  bar  fixed  at  the  outer 
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wall  ol  the  test  box.  In  other  experiments  with  a  window  of  5  cm 
thickness,  the  plug  was  allowed  to  move  out  of  the  hole.  In  this 
case,  the  high-speed  movie  camera  may  look  at  the  region  Just 
below  the  charge,  which  was  covered  by  a  lock  bar  in  previously 
mentioned  experiments.  To  be  sure  that  the  plug  does  not  make  any 
damage,  a  rope  kept  the  plug  a  on  a  short  rotational  trajectory 
to  the  floor - 

Small  plastic  discs  and  colored  sand  were  used  to  mark  the 
subsurface  deformation.  Plastic  discs  were  used  to  measure  final 
ejecta  velocities,  colored  sand  quadrangles  were  used  to  measure 
trajectories  of  corner  points.  Pig.  5.1  illustrates  the  motion  of 
4  near-surface  plastic  markers  in  a  quarter-space  experiment - 

The  first  experiments  were  conducted  with  a  modestly  dense 
wet  sand.  Experiments  shown  two  stages  of  the  cavity  growth:  At 
the  first  stage  the  crater  grown  as  a  self-similar  hemisphere. 

The  apparent  center  of  the  hemisphere  was  at  the  pre-event 
surface  are  a  little  bit  above.  With  this  geometry  the  cavity 
grown  up  to  the  maximum  depth.  At  the  second  stage,  the  bottom  of 
the  cavity  rest  at  the  place,  but  side  parts  of  the  target 
continued  to  move,  creating  so-called  "wings"  of  a  crater.  The 
final  depth-diameter  ratio  was  close  to  0.5  for  wet  sand.  Fig. 

5.2  illustrates  the  observed  crater  growth. 

Trying  to  describe  the  motion  around  the  cavity  during  the 
first  stage,  we  went  to  the  conclusion  that  the  radial  velocity 
at  all  cavity  surface  may  be  supposed  as  a  constant  value 
decreasing  with  time. 

In  parallel  to  these  experiments,  we  analyzed  the  flow  of 
water  due  to  surface  explosion.  At  the  late  stage,  when  the 
incompressible  model  may  be  used,  the  cavity  shape  in  water  was 
very  close  to  a  hemisphere,  what  was  very  instructive  being 
compared  with  surface  bursts  in  sand  targets.  We  constructed  the 
self-similar  potential  flow,  using  the  series  of  harmonic 
functions,  revealed  close— to— power  law  spatial  decay  of  flow 
velocities. 

Combining  these  two  data  sets,  we  got  two  main  supposition 
for  Z-model.  If  one  uses  the  polar  coordinates  measuring  radial 
distances  from  the  center  and  the  angle  from  the  down-ward 
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vertical  line,  these  suppositions  may  described  as  follows:!. 
Radial  velocity  is  a  power  function  only  of  radial  distance. 2. 
Tangential  velocity  is  defined  by  integral  balance  equations  (lor 
simplest  case  -  the  balance  of  mass). 

Por  target  materials  of  constant  density  ( "incompressible 
media")  it  is  easy  to  derive  the  velocity  field  from  above 
mentioned  conditions. 

We  use  the  polar  coordinate  system  R,  iff,  tp.  The  5  =  0  axis 
presents  the  downward  vertical  direction,  R  is  the  distance  from 
the  coordinate  center,  tp  is  an  angle  at  horizontal  plains.  As  we 
treat  axis-symmetrical  problems  here,  all  values  do  not  depend  on 
p.  Por  such  a  coordinate  system  depth,  2,  is  expressed  as 

z  =  R  cos  5 

and  radial  distance  from  the  vertical  axis,  r,  is  expressed  as 

r  =■  R  sin  ff 

After  the  explosion  deposition  of  energy  at  the  point  R=0,  shock 

waves  go  out  from  the  explosion  source.  The  target  material 

behind  the  wave  front  begins  to  move.  When  the  shock  waves  go  far 

from  the  center  and  decay,  the  zone  of  Inertial  motion  is  formed 

at  the  excavation  area.  The  velocity  field  in  this  area  may  be 

approximately  described  as  follows. 

At  any  point  the  velocity  vector  may  be  presented  as  the 

vector  sum  of  radial,  and  tangential,  components. 

According  to  Z-model,  for  near-surface  explosions  the  value  of  v„ 

H 

does  not  depend  on  the  angle  ff  and  may  be  described  as  a  power 
function  of  R: 

(5.1) 

The  condition  of  incompressibility  in  spherical  coordinates  is 
written  as: 

^  -  -4-  *  3vn  4^  <■"8  °  <5-a) 

Substitution  of  (5*1)  to  (5-2)  gives  the  equation  for  the 
tangential  velocity  component,  which  has  a  solution: 

=  (Z-2)  (1-  cos  ^)/stn  5  (5.3) 

So  from  assumption  (5-1)  and  (5*2)  we  get  the  following  result: 

The  tangential  velocity  component,  depends  on  radial  distance 
in  the  same  way  as  Vj^t  and  at  a  given  R  gradually  increases  from 
zero  value  at  the  axis  of  symmetry  5=0  to  the  maximum  value  at 
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the  free  surface  ^=90°  (If  the  coordinate  center  is  situated  at 
the  free  surface ) . 

Por  the  flow  field  expressed  by  (5.1)  and  (5-3)  it  easy  to 
get  trajectories  of  particles: 

/?  =  [(1  -  cos  9)/(1-cos  (5.4) 

where  and  are  coordinates  of  a  "pin”  point  -  a  point  which 
fixes  a  given  trajectory. 

The  Eq.  (5.4)  clearly  demonstrate  the  the  shape  of 
trajectories  does  not  depend  on  time,  so  trajectories  are 
stationary  streamlines.  Note  that  he  motion  of  particles  along 
these  trajectories  may  be  non-stationary  one-  Rotating  of  two 
streamlines  around  the  vertical  axis  of  symmetry  gives  the 
so-called  stream  tubej  which  define  the  excavation  motion  of 
target  material. 

Por  some  problems  it  is  suitable  to  choose  as  the  "pin” 
point  the  point  of  a  trajectory  and  the  initial  target  surface 
level  intersection:  R^  “  ^e’ 

According  to  Z-model  the  motion  of  material  ejected  from  a 
crater  is  described  as  motion  along  streamlines,  defined  with  the 
equation 

R  =  “Tq  (1  -  cos  (5.5) 

Here  r^  is  the  radial  distance  where  a  streamline  crosses  the 

initial  target  surface  level  (ITSL),  the  coefficient  Z  is  in  the 

range  from  2.5  to  3  for  most  cases.  A  pair  of  two  neighbor 

streamlines  which  cross  the  ITSL  at  radial  distance  r  and  r  +dr 

e  e  e 

define  a  stream  tube  with  volume 

'■  2-2 

^  ”  [—  J  <5.6) 

The  total  ejecta  volume  from  a  crater  with  the  radius  R^  at  the 
ITSL  is 

^  2-2 

V  =  (2/3)  n  -  flf  (5.7) 

I  Z+1  J  ^ 

Eqs.  (5.6)  and  (5-7)  are  derived  for  the  model  with  the  center 
(/?=0,  5=0)  at  the  ITSL.  The  similar  equations  for  model  with 
burled  center  may  be  found  elsewhere  (sf.  Croft,  1981). 

The  ejecta  velocity  for  gravity  craters  may  be  estimated 
from  rather  simple  consideration  (Ivanov,  1979)-  The  potential 
energy  of  the  streamline  in  the  gravity  field  with  the  free-fall 
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acceleration  g  is 


dP^  = 


P  g  d7^ 


CZ-2)(Z+1) 


2  Z  (Z+2) 


(5-8) 


If  one  approximate  the  ejecta  initial  velocity,  as  an 
exponential  function  of 

“e  -  "r  <5-9) 

where  is  the  coefficient  which  may  be  called  as  "the  ejecta 
velocity  at  the  crater  edge".  Por  the  simplest  estimation  m  may 
be  treated  as  a  constant  in  the  range  from  2.5  to  3. 

The  kinetic  energy  at  the  moment  when  all  the  material 
ejected  along  the  given  stream  tube  is  raised  above  the  ITSL  may 
be  estimated  as 


c3/r  -  p 


ruf/2) 


The  ejection  process  is  arrested  at  r  =  where 

cK-dP 


or 


P 


p  (vy2)  ^  p  g 


r2-2)(Z+1) 


2  Z  (Z+2) 

The  value  of  may  be  expressed  from  the  Eq.  (5-12)  as 


where 


k  = 


k  g  R^ 

''  CZ-2)(Z+1) 


(5.10) 


(5.11) 


(5.12) 


(5.13) 


Z  (Z+2) 


e’ 


(5-14) 
as  a  function 


Eqs.  (5.9,  13,  14)  define  the  ejecta  velocity,  v 
of  the  radial  distance,  r^,  of  the  point  of  ejection  trough  the 
ITSL.  The  angle  of  ejection  at  the  ITSL  is  depended  on  Z  value  as 

<p  =  tan"”'  iZ-2)  (5.15) 

The  ejecta  motion  above  the  ITSL  may  be  modeled  with  the 
well-known  equation  of  a  a  free  body  flight  in  the  gravity  field. 
The  horizontal  component  of  the  velocity  is  constant  up  to  the 
moment  of  deposition  beyond  the  crater  rim  and  is  equal  to 


V 


eh 


Ug  sin  <p 


(5.16) 


The  horizontal  velocity  of  the  main  mass  of  ejecta,  deposited  at 
the  rim  may  be  estimated  as  sin  (p.  Further  from  the  rim  the 
horizontal  velocity  gradually  increases. 

Beyond  the  crater  edge  (for  >  R^)  the  material  in  stream 
tubes  does  not  ejected  but  displaced  into  a  new  position.  The 
part  of  the  material  raised  above  the  ITSL  constructs  the 
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so-called  structural  uplltt  at  the  base  of  a  crater  rim-  The  real 
motion  in  this  region  is  rather  complicated,  but  at  the  first 
approximation  it  may  be  described  in  terms  of  Z-model  as  follow. 

Supposing  the  constant  density  of  the  target  material  one 
may  derive  that 


(' 


(5-17) 

whera  is  the  radial  distance  to'  the  given  streamline  at  the 

ITSL,  is  the  final  radius  of  the  crater  profile  point  of 

intersection  with  the  given  streamline,  Rj,  is  the  radial 
coordinate  of  top  point  of  the  streamline  at  the  uplifted  target 
surface.  Lower  trajectory  at  the  Fig.  5-3  illustrates  the  model 
treatment  of  under-crater  mass  motion. 

Using  the  geometric  relationship  (5-5)  one  may  calculate  the 
cartesian  coordinates  of  the  uplifted  target  surface.  At  the 
boundary  between  ejected  and  uplifted  material  (r  =  a.  =R^.) 

R,  ^  '(5.18) 

And  the  Cartesian  coordinates  of  the  uplifted  rim  crest  are 

^ur 


R^  sin  «  ^ 
t  ur 


vr 

cos  9. 


Rt  cos 


(5.19) 


.2-2 


=  r  -  CRj^R^r 

Eqs.  (5.19)  give  the  maximum  estimate  of  the  uplifted  rim  height, 
probably  never  realized  in  nature  due  to  complex  material  flow  at 
the  crater  rlra  and  immediate  subside  of  the  rim  for  complex 
craters.  Nevertheless  Eqs.  (5.18)  and  (5-19)  allow  to  make 
first-order  estimates  of  the  uplifted  rim  profile. 

For  Z=3  one  can  derive: 


(jt/6) 


.3  . 


t  ’ 

=  (n/2)  rdr 


e* 


V 


R 


=  /  (4/15)  g 


R. 


t’ 


ft  =  4/15 


<P 


45°; 


ur 


(2/15)  g  R^ 

=  1.168  R^ 


=  .365  / 


g  Ri 


(5.20) 


where  is  the  total  ejected  volume  from  the  crater  of  an 
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apparent  radius  (37^  is  the  volume  ejected  through  the  stream 
tube  ot  dTg  width  at  the  surface,  tp  is  the  ejection  angle,  is 
the  characteristic  ejection  velocity  near  the  gravity  crater 
edge,  is  the  horizontal  component  of  the  ejecta  velocity,  g 
is  the  gravity  acceleration,  the  radius  of  the  structural 
Uplift  rim  crest,  z^  is  the  amplitude  of  the  structural  uplift 
of  the  pre-event  surface.  These  values  are  widelyused  in  various 
applications . 


5.3.  Using  of  the  Z-model 

During  70th  and  80th  we  used  the  model  to  estimate  various 
application  of  surface  burst  cratering.  As  some  of  results  has 
been  published,  and  the  same  approaches  has  been  found  in  the 
American  report  [8],  here  we  shortly  list  several  main 
applications: 

1.  Cratering  in  the  media  with  an  ideal  (von  Mises)  plasticity 
was  studied  with  the  Z-model  in  [11].  Analogous  study  has  been 
found  in  [8]. 

2.  The  Z-model  helped  us  to  introduce  the  separation  of  gravity 
and  strength  regimes  of  cratering  [9,  12,13].  This  approach  was 
Independently  elaborated  by  R. Shmidt  and  K.Holsapple  [14,  15]. 

3.  Using  Z-raodel  we  elaboration  a  functional  form  to  construct  a 
form  of  a  wide-range  scaling  law,  recently  used  for  planetary 
applications  [16]. 

4.  The  mass  balance  based  on  Z-model  helped  us  to  create  the 
model  of  the  ejecta  layer  deposition  [9].  It  become  obvious  from 
tills  study  that  rim  profiles  for  all  gravity  craters  should  to  be 
geometrically  similar.  The  same  result  was  reproduced  by  Housen 
o.a.  [17]. 

As  these  results  were  published  (sometime  in  a  short  form), 
it  seems  practical  to  skip  them  here.  We  can  return  to  these 
questions  if  anybody  will  be  Interested  to  discuss  them.  In  this 
report  we  concentrate  at  the  problem  of  Z-raodel  with  dilatancy, 
which  analog  we  never  saw  at  the  literature.  The  short 
description  of  the  model  was  published  in  [18]. 
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5.4.  Z-model  with  shear  bulking. 


More  than  100  years  ago  Reinolds  noted  that  shear 
deformation  of  grain  media  is  accompanied  with  the  increase  of 
specific  volume.  This  phenomenon  was  named  as  bulking  or 
dllataney.  Recently  the  dilatancy  was  under  investigation  in 
connection  with  stress  wave  and  general  rock  mechanics  [20].  In 
this  section  we  demonstrates  how  the  Z-model  may  be  extended  to 
describe  a  cratering  flow  with  dilatancy,  i.e.  when  the  volume  of 
porous  space  in  the  target  media  may  increases  due  to  shear 
strain.  To  do  it  we  need  to  sacrifice  with  a  generalization  of 
Z-model  for  arbitrary  crater  shape.  The  value  of  void  volume  may 
characterize  to  some  extent  the  level  of  damage  of  rocks  under 
the  crater. 

The  simplest  representation  of  dilatancy  is  the  kinematic 
relation,  described  proportionality  of  the  volume  deformation 

rate,  I  =  Y/V  =  v  v  (V  -  is  specific  volume  of  material),  to  the 
second  invariant  of  the  strain  rate  tensor. 


where  -0,5  (dv^/dXj  +  dVj/dx^). 
Por  such  a  media  the  mass  conservative  law  is  expressed  as: 
1  a  P  1 

_  Vj.)  +  — 

jyC  r> 
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dv 


a  {  0,5  [( 
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R  Gin  5  35 
1  31/ „  V 
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V, 


/  n  p  /  ^V n  ^ 

+  -  ctg  d  -  ^  ^  )  J  +  3  (-^  ^ 


(5.21) 


where  d  is  the  coefficient  named  as  dilatancy  coefficient . For  the 
power  law  decay  of  radial  velocity  with  a  distance  (Eq.  5.1),  the 
above  equation  is  simplified  to  the  form  below: 

-(^'-2)  +  /'  +f  ctg  d  = 

-  d  (0,5  [(Z+1+f)^  +  (Z+Uf  ctg  5;^  +  (f  -  f  ctgd)^]  + 
i  0,75  (Z+1)^  (5.22), 

where  function  /(5)  describe  the  angle  dependence  of  a  tangential 
velocity: 
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•  fC^)t  and  /'=  df/dd. 

Itajectorles  of  particle  motion  now  is  described  as 

^  -r 

R/R.  ^  exp  C  J  f  ’  d^) 

'  a 


(5.23) 


which  replace  the  Eq.  5-5- 

_One  can  see  that  for  hemispheric  cavities  the  including  of 
dilatancy  in  the  linear  approximation  conserves  the  main  result 
of  the  traditional  Z-model :  any  particle  moves  along  a  trajectory • 
(a  streamline )  which  does  not  change  in  time .  Nevertheless  the 


specific  volume  of  any  element  INCREASES  due  to  porous  space 
increases  by  dilatancy: 


y/7  =  V  u  =  (A(t)/R)  [-(Z-2)  +  f'  +  f  ctg  (5.2A) 

where  A(t)  -  is  a  function  of  time,  depended  on  the  cavity  growth 
rate. 

Using  of  an  exact  solution  for  the  angle  function  /(3)  from 
the  Eq.  5-22  is  inconvenient  as  it  destroy  the  main  advantage  of 
Z-model  philosophy  of  simple  analytical  estimates.  However  our 
analysis  shown  that  the  f unction/r^ ^  may  be  affordably 
approximated  in  the  familiar  form  of  the  Eq.  5.5  if  one  use  the 
approach  described  below.  One  can  see  that  the  main  effect  of 
dilatancy  for  a  given  Z  value  is  the  more  steep  curyation  of 
streamlines  to  the  surface  in  comparison  with  the  zero  dilatation 
((3=0).  As  a  rule  in  experiments  one  can  measure  deformation  field 
much  more  easier  than  a  velocity  field.  Then  one  can  split  the 
single  value  of  Z  in  the  traditional  model  into  2  values:  and 

Zy,  where  Z^  will  describe  trajectories  in  the  form  of  Eqs.  5.5, 
while  Zy  value  will  replace  the  "old”  Z-value  in  the  velocity 
field  description  in  Eqs.  5-1  and  5.3.  The  relation  ship  between 
Z^  and  Zy  may  be  found  from  the  following  consideration. 

Let  the  stream  tube  has  a  form  of  Eq.  5.5  with  the  parameter 

Z=Z^: 

R  ^  Rj  [(1  -  cos  9)/(1-C03  (5.25) 

and  the  velocity  field  is  described  by  Eqs.  5.1,  5.2  with  the 
parameter  To  fit  the  streamline  one  needs  to  determine  the 

value  of  f'(9-0).  Using  the  limit  transition  one  can  find 
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that/’r^=0;=- 


(Zy-2)  +  do(Z^-hl) 

'  5=a 


Por  zero  dilatancy,  £3=0  (as  in  the  standard  Z-raodel 

f^(d=0)  =  Z^-2 

The  initial  derivative  of  f(^)  functions  coincide  if 

2o(d+1> 

2-d 


(5.26) 


For  simple  radial  flow  for  buried  explosion  Z^  =  2  and  the  Eq. 
5.26  has  a  form 


2-d 

“  T~nj 


(5.27) 


Previously  the  same  relationship  was  derived  in  [  20  ] ,  and  was 
successfully  used  for  the  buried  explosion  analysis. 

As  one  can  find  better  way  to  choose  proper  pairs  of  Z^  h 
the  Eq.  5.27  appears  to  be  enough  correct  for  practical 
estimates. 

The  first  check  of  the  model  we  tried  to  do  with  the  data  by 
Piekutowski  [19].  who  measured  the  dilatancy  of  dry  sand  under 
the  crater.  Pig.  5.4  a  and  b  compare  the  volume  strain 
distribution  under  the  crater.  The  best  fit  occurs  Just  under  the 
crater.  The  model  does  not  predict  enhanced  near-surface  bulking, 
possibly  due  to  air  blast  or  some  spallation  effects  of  stress 
waves.  Anyway  the  fit  may  be  recognized  as  relatively  good, 
taking  into  account  the  simplicity  of  the  model  we  used.  Original 
data  similar  to  published  in  [19]  may  give  the  new  challenge  to 
test  and  to  refine  the  model. 


5.5.  An  example  of  mass  balance  for  surface  nuclear  events. 
Uelow  we  present  an  example  of  nuclear  near  surface  burst 
cratering  volumetric  analysis  for  event  #4  and  #11  listed  in  the 
Section  3  of  a  current  report. 

Integration  of  the  ejecta  deposits  thickness  around  craters, 
31sted  in  the  Section  3,  allows  to  estimate  the  low  limit  of  the 
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total  volume  of  ejecta.  The  fitting  of  these  data  with  the 
Z-model  demonstrates  the  typical  value  of  Z^=2.5±0.1.  Due  to  the 
lack  of  direct  measurements  during  surface  events,  it  is  hard  to 
find  the  way  to  check  these  values  with  any  direct  transient 
displacement  data.  Only  profiles  of  a  structural  uplift  allows  to 
check  these  values  by  restoring  the  uplift  profile  from  the 
Z-raodel  streamline  geometry  and  observed  shape  of  a  crater.  Por 
current  estimates  we  supposed  the  center  of  coordinates  to  be  at 
the  level  of  the  target  surface.  It  is  possible  in  principle  to 
vary  the  position  of  the  coordinate  center  for  better  fitting. 

To  define  the  crater  profiles,  we  used  a  scanner  to  digitize 
crater  profiles,  shown  on  pictures  in  the  Section  3.  The  data  on 
direct  measurements,  also  displayed  in  the  Section  3,  were  used 
to  control  the  scale  of  digitized  pictures.  Por  any  point  at  the 
crater  profile,  the  streamline  was  was  traced  above  the  prelrapact 
process.  The  mass  conservation  (with  a  proper  volume  increase) 
determined  the  final  position  of  the  structural  uplift.  Closer  to 
the  crater  rim,  the  theoretical  rim  profile  is  too  steep  to  be 
stable,  so  partial  down  slide  seems  to  be  unavoidable. 

It  is  become  obvious  very  soon  that  the  value  of  Z 
parameters,  estimated  from  total  ejecta  volumes  are  too  low  to 
fit  the  structural  uplift  and  the  observed  crater  profile.  At  the 
same  time,  the  mean  density  of  the  uplift  is  definitely  smaller 
than  initial  target  material  density.  So  the  supposition  of 
dilatancy  of  a  matter  during  cratering  becomes  to  be  natural. 
Unfortunately,  the  theoir  described  at  the  previous  section  works 
exactly  only  for  hemispheric  cavities.  If  the  shape  of  the  cavity 
differs  from  the  hemisphere,  it  means  that  shear  strain  may  be 
larger  than  for  a  hemispheric  cavity  (neighbor  streamlines  have 
larger  relative  displacement ) . 

Nevertheless  for  simple  estimates  we  used  the  model 
relations  Eqs.  5*4  (with  Z=Z^)  and  5-26  to  describe  the  cratering 
flow  with  dilatancy.  In  such  an  approximation  the  modest  value  of 
the  dilatancy  coefficient  d  ss  0.2  produce  much  better  fitting 
than  the  standard  Z-raodel.  As  uplift  fitting  allow  to  estimate 
the  "trajectory"  value,  the  "velocity"  Z-value,  Z^,  is 
calculated  from  Eq.  5*26. 
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Pigs-  5-5  and  5-6  demonstrates  the  best  fitting  we  found  for 
4  cross  section  of  each  of  craters. 

Por  events  #4  and  #11  the  best  value  of  was  found  as  2.6 
and  2.6?  correspondingly  (which  coincides  with  the  total  ejecta 
volume  estimates).  Por  d=0.2,  the  value  of  Zy  is  about  2.0  —  much 
lower  than  traditional  fitting  of  crater  computational  velocity 
field!  2.7  to  3. 

The  dilatancy  is  naturally  suppressed  by  high  pressures,  as 
well  as  the  dilatancy  stop  to  be  linear  one  as  the  material  begin 
to  be  very  loose.  So  the  range  of  pressure  and  volume  strain 
exist  where  the  dilafancy  really  can  play  some  role.  The  effect 
was  thoroughly  Investigated  for  stress  wave  decay ,  but  it  is  very 
Important  to  find  the  proper  place  of  the  dilatancy  in  cratering 
flow  field  description. 
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Fig.  5.1.  a).  An  example  of  a  quarter-space  experiment.  The  high 
speed  move  frames  was  overlapped  to  plot  ejecta  trajectories. 
Plastic  disks  embedded  to  the  wet  sand  are  ejected  with  decreasing 
with  the  ejection  distance  velocities,  (b)  Z-model  representation 
of  the  same  experiment.  The  material  ejected  along  a  stream  tube  a 


PiekuloTTski,  1900 


Fig. 5.-!.  ia)  Experimental  data  [19]  illustrating  the  sand. 
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Fig.  5.5.  Structural  uplift  for  the  surface  burst  =4  in 
comparison  with  Z-model  estimates. 


Fig.  5.6.  Structural  uplift  for  the  surface  burst  #11  in 
comparison  with  Z-raodel  estimates. 


Held  data  SU 
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Appendix 

Terminology  and  Nomenclature 

Below  we  explain  some  detailes  of  terminology  for  sedimentary  rocks  used  in  the 
presented  report. 

a.ieBpiiT_(alevrite)  -  pbixnas  ocaaoMHas  nopoaa,  no  cocrasy  npoMexjTOMHaa  Me2ca>’ 
necMaHbiNm  h  rmiHUCTbiNOi  noponaMw  c  pcuiMepoM  rnaBHofi 
Maccbi  aepcH  0.01  -0.1  mm 

aaeBpoaiiT  (alevrolite)  -  cueMCHTiipoBaHHbifi  aneBpitr;  cucMCHTHpoBaHHbiH 
TOHKojrncnepcHbiH  Marepnan  (6o.iee  50  %  Macruu  0.01-0.  I.mm 
FLioTHOCTb  2.68-2.75  p/c.mS, 
npoHHOCTb  Ha  cxaTHe  49-78  Mil  a 
npoMHOCTb  Ha  pacTsxeHiie  5.9-8. 8  MHa 
CKOpOCTb  npOflOJIbHbIX  BonH  4500-5500  .m/c 

apnmnirr  (argilite)-  ocaao'-maa  ropnaa  nopoaa,  oOpaaoBasmaaca  b  peayabTaTe 

yruTOTHeHHa,  oOesBoacHBawia  ii  ueMeHiaunH  ituih;  ornHMacTca  or  rntm 
6oabmeH  TBepaocxbio  ii  HecnocoGHocibio  pa3.MOKaTb  b  Boae 

c)Tiecb  (sandy  loam,  sandy  soil)-  necMHO-nuiHiicTaa  nopoaa,  coaepacamaa  10-20% 
npfLVieCH  mHHHCTbIX  HaCTHU  cjTiecb, 
cjnecManaa  noMsa  =  sandy  clay 

c\T.iii]iOK,  cyTaHHucraa  noHsa  =  sandy  loam 

cjT.iiiHOK  -  (loam,  loam  soil)  rawHa  hwskoh  nnacTHMHOCTw  c  npuMccbio  30-40% 
necKa 

rmiHiicTbin,  cjTraiKiicTbin,  McpreabHbin  =  loamy 
6^75111  cyraiiHOK  =  silt  loam  • 

raiiHa  =  shale 

sandy  shale 
carbonaceous  .shale 


r.iHHiicTbin  caaneu,  caanuesaa 
necHaHo-raiiHiicTbiH  caaHeu  = 
jT-iiicTbiu  caaneu  = 


.  Grain  size: 


Psepbites  (rudites,  gravel,  rubble,  boulder,  pebble)  are  clusiic  sedimentswith  more  than 
50%  (by  weigth)  of  the  constituents  >  2  mm  in  diameter. 

Psammites  (arenites;  sand)  have  more  than  50%  of  their  grains  within  the  sand  fraction 
(2  to  0.002  mm) 

Pelltes  (argillites;  clay,  mud)  are  restricted  to  grain  size  <  0.002  mm.  A  transitioana! 
grain  size  fraction  within  the  psammite  fraction  is  that  of  silt  (from  0.06  to  0.002  nun) 


Approximate  correspondence  of  English  and  Russian  nomenclature  of  sedimentary' 
rocks:  English  nomenclature  according  to  [3];  Russian  nomenclature  (I):  according  to 
{!];  (II):  according  to  geologic  description  in  primary  reports. 


Particle 

Diameter, 

mm 

<1/256 
(<  0.004) 

1/256  to 

1/16 

(0.004  to 
0.06) 

1/16  to  2 
(0.06  to  2) 

2  to  4 

4  to  64 

64  to  256 

>256 

Particles 

Clay 

AJevrite 

(SUt) 

Sand 

Gra\'el 

Pebble 

Cobble 

Boulder 

Kind  of 

rock 

Clay 

Argilite 

AJevrite 

Alevrolite, 

Siltstone 

Sand, 

sandstone 

Conglomerate,  Breccia 

MaCTHUN 
MM  (I) 

AnespHT 
0.01  to  0.1 

IlecoK 

0.1  to  1 

rpaBtm 

1  to  10 

rajibKa 

10  to  100 

sanyHbi 

100  mm  to 

10  m 

HacTHUbi 
MM  (IT) 

apecBa 

2  to  20 

uieOeHb 

20  to  200 

BanyHbi 

>200 

flopoabi 

rniiHa 

Apirumnr 

AncBpirr, 

aneBpomiT 

( 

riecoK 

necMaHtiK 

KoHrno.MepaT,  SpcKmis 

apecBa  2-20  mm  =  pebble 
ineOcHb  20-200  mm  =  rubble 

Both  may  be  named  as  breccia  with  different  fragment  size. 
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DEPARTMENT  OF  THE  AIR  FORCE 

AF  TECHNICAL  APPLICATIONS  CENTER 
ATTN:  ROBERT  BLANDFORD 

AFSOR/NP 

ATTN:  COL  JERRY  PERRIZO 

AIR  FORCE  TECHNICAL  APPLICATIONS  CTR 
ATTN:  FRANK  PILOTTETT 
ATTN:  GEORGE  ROTHETTR 

AIR  UNIVERSITY  LIBRARY 
ATTN:  AUL-LSE 

GEOPHYSICAL  LABORATORY 

ATTN:  JAMES  LEWKOWICZ 
ATTN:  XO 

HQ  497 IG/INOT 
ATTN:  INT 

DEPARTMENT  OF  ENERGY 

DEPARTMENT  OF  ENERGY 
NEVADA  OPERATIONS  OFFICE 
ATTN:  OTIS  D  H  MARTIN 

DEPARTMENT  OF  ENERGY 
ATTN:  DR  CV  CHESTER 
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LAWRENCE  LIVERMORE  NATIONAL  UB 
ATTN:  RDONG 
ATTN:  DB  CLARKE 
ATTN:  WILLARD  J  HANNON  JR 
ATTN:  G  POMYKAL 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  BOB  DEUPREE 
ATTN:  M  HENDERSON 
ATTN:  PAUL  JOHNSON 
ATTN:  FREDAPP 
ATTN:  TKUNKLE 
ATTN:  R  P  WEAVER 
ATTN:  DSTROTTMAN 
ATTN:  REPORT  LIBRARY 
ATTN:  S  TAYLOR 
ATTN:  R  KIRBY 
ATTN:  TECH  LIBRARY 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  J  M  MCGLAUN  6444 
ATTN:  MIKE  FURNISH 
ATTN:  T  K  BERGSTRESSER  DIV  931 1 
ATTN:  TECH  LIB  3141 
ATTN:  9312  CW  SMITH 
ATTN:  9312  D  GARBIN 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  KATIE  POLEY 
ATTN:  OSWR/NEDL  TURNBULL 
ATTN:  OSWR/NED  5S09  NHB 

U  S  ARMS  CONTROL  &  DISARMAMENT  AGCY 
ATTN:  A  LIEBERMAN 

US  GEOLOGICAL  SURVEY  NATIONAL  CENTER 
ATTN:  W  LEITH 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORP 

ATTN:  LIBRARY  ACQUISITION 

APPLIED  RESEARCH  ASSOCIATES,  INC 
ATTN:  CJ  HIGGINS 

APPLIED  RESEARCH  ASSOCIATES,  INC 
2CYATTN:  SBLOUIN 

APPLIED  RESEARCH  ASSOCIATES,  INC 
ATTN:  R  FRANK 

BDM  FEDERAL  INC 

ATTN:  EDORCHAK 
ATTN:  J  STOCKTON 

BOEING  TECHNICAL  &  MANAGEMENT  SVCS,  INC 
ATTN:  ROBERT  M  SCHMIDT 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 
ATTN:  T  AHRENS 

CARPENTER  RESEARCH  CORP 
ATTN:  HJ  CARPENTER 


ENSCO,  INC 

ATTN:  DOUGUSBAUMGARDT  - 

HORIZONS  TECHNOLOGY,  INC 
ATTN:  BLEE 

IIT  RESEARCH  INSTITUTE 

ATTN:  DOCUMENTS  LIBRARY 
ATTN:  M  JOHNSON 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  CLASSIFIED  LIBRARY 

JAYCOR 

ATTN:  CYRUS  P  KNOWLES 

KAMAN  SCIENCES  CORP 
ATTN:  DCAYNE 
ATTN:  LIBRARY 

KAMAN  SCIENCES  CORP 
ATTN:  DASIAC 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 

KTECH  CORP 

ATTN:  F  DAVIS 

LACHEL  AND  ASSOCIATES,  INC 
ATTN:  J  BECK 

LOGICON  R  &  D  ASSOCIATES 
ATTN:  C  K  B  LEE 
ATTN:  D  SIMONS 
ATTN:  LIBRARY 

LOGICON  R  &  D  ASSOCIATES 
ATTN:  J  RENICK 
ATTN:  L  GERMAIN 

LOGICON  R  &  D  ASSOCIATES 
ATTN:  D  PIEPENBURG 

MAXWELL  LABORATORIES  INC 
ATTN:  DR  JL  STEVENS 
ATTN:  DR  K  L  MCLAUGHLIN 
ATTN:  DRCDISMUKES 
ATTN:  KDPYATTJR 
ATTN:  RUFRENZ 
ATTN:  S  PEYTON 

MAXWELL  LABORATORIES,  INC 
ATTN:  J  MURPHY 

MCDONNELL  DOUGLAS  CORPORATION 
ATTN:  R  HALPRIN 

MISSION  RESEARCH  CORP 
ATTN:  TECH  LIBRARY 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  H  BRODE 

S-CUBED 

ATTN:  JACK  MURPHY 


Dist-2 
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SCIENCE  APPLICATIONS  INTERNATIONAL  CORP 
ATTN:  RSCHLAUG 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  G  EGGUM 
ATTN:  G  PHILLIPS 
ATTN:  H  WILSON 
ATTN:  L  SCOTT 
ATTN:  M  MCKAY 
ATTN:  TCBACHEJR 
ATTN:  TECHNICAL  REPORT  SYSTEM 
ATTN:  WWOOLSON 


TITAN  CORPORATION  (THE) 

ATTN:  LIBRARY 
ATTN:  S  SCHUSTER 

TRW  INC 

ATTN:  TIC 

TRWS.l.G. 

ATTN:  NORMAN  LIPNER 

TRW  SPACE  &  DEFENSE  SECTOR 
ATTN:  W  WAMPLER 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  COCKAYNE 
ATTN:  WLAYSON 

SOUTHERN  METHODIST  UNIV 
ATTN:  DR  BRIAN  STUMP 

SRI  INTERNATIONAL 

ATTN:  A  FLORENCE 
ATTN:  DR  JIM  GRAN 

TECH  REPS.  INC 

ATTN:  FMCMULLAN 


WASHINGTON,  UNIVERSITY  OF 
ATTN:  J  KATZ 

WEIDLINGER  ASSOC,  INC 
ATTN:  H  LEVINE 

WEIDLINGER  ASSOCIATES,  INC 
ATTN:  TDEEVY 

WEIDLINGER  ASSOCIATES,  INC 
ATTN:  I  SANDLER 
ATTN:  M  BARON 

FOREIGN 


TERRA  TEK,  INC 

ATTN:  S  GREENE 

TITAN  CORPORATION 

ATTN:  J  THOMSEN 


INSTITUTE  FOR  DYNAMICS  OF 
aCYATTN:  B  IVANOV 
2CYATTN:  LPERNIK 
2CYATTN:  VADUSHKIM 


